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1.   INTROKJCTIQK 

Every  industrial  chenical  process  is  desisned  to  produce  econoraically  a 
desired  product  fl-om  a  variety  of  starting  materials  throuch  a  succession  of 
treatment  steps.  "Hie  raw  materials  underco  a  number  of  physical  treatment 
steps  to  put  them  in  a  form  in  ^diich  thecir  can  be  reacted  chemically.  They 
then  pass  th.rough  the  reactor.  The  products  of  the  reaction  must  then  undei^o 
further  physical  treatment  (separations,  purifications,  etc.)  for  the  final 
desired  product  to  be  obtained. 

The  chemical  treatment  step  is  oftm  the  heart  of  the  process,  the  thing 
that  makes  or  breaks  the  process  economically.  Design  of  the  reactor  is  no 
routine  matter,  and  a  great  variety  of  designs  can  often  be  proposed  for  a 
process.  In  searching  for  the  optimum  design  it  is  not  just  the  cost  of  the 
reactor  that  must  be  minimized. 

Most  chemical  plants  consist  of  a  large  number  of  inter-connected  process, 
units  x^ich  together  fora  a  topologically  conplex  system  (a  systeci  with 
branches  and  loops).  The  chemical  reactors  are  often  arranged  in  such  a  way. 
To  obtain  -toe  optimal  design  of  a  geometrically  coji^Jlea:  process,  we  have  to 
select  the  optimal  process  parameters,  and  operating  variables  according  to  a 
certain  performance  criterion  (objective  fimction).  Ilany  r^thods,  e.g.  differ- 
ential calculus,  calc'jlus  of  variations,  linear  and  non-linear  prc^ramming, 
dynamic  programming,  and  the  maximum  principle  have  been  developed  to  solve 
these  optimization  problems.  Most  of  them,  hotwver,  are  not  suitable  for 
practical  use  ^en  the  processes  become  geometrically  ccMplex  and  the  pirocesa 
variables  increase,  'dynamic  programming  and  the  maximum  principle  were  derived 
and  have  been  greatly  extended  to  attack  such  gecaastricaUy  complex  and  miltl- 
variable  processes. 


Ihe  raaxiimm  principle  was  first  developed  for  slraple  contimaous  processes 
by  the  Russian  mthematician  Pantr/agin  j  Ij.  Katz  ^Z,   3J  J^as  presented  a 
discrote  version  of  the  raaxiimm  priJiciple  for  simple  hcasK^oneous  gystens. 
Jackson  j  ^i  j  has  pointed  out  that  Kata's  algorithm  is  only  true  as  a  "weak 
algorithEi".  This  so-called  "local  weak  maxiraun  principle"  <Mr  "stationary 
inrinciple"  on!;?'  relates  stationaiy  values  for  the  discrete  case.  The  exten- 
sion of  Katz's  results  by  Fan  and  V&ng  \5,   6,  7 J  to  gystens  topologicaTly  raore 
c<aJiplex  than  a  sij^lo  sequential  chain  is  also  onl/  a  vmsk  alcorithra  relating 
stationary  values.  Fan  and  l-feng  j  8  j  have  Qcctendod  and  generalized  the  ccai- 
tinuoas  raaxiimxm  principle  so  that  it  can  be  used  to  optiraize  topologically  or 
geosnetricallc^''  ootaplex  continuous  processes.  Several  other  successful  att«apts 
to  "generalize"  the  discrete,  the  continuous,  and  the  composite  versions  of 
the  maximum  principle  for  ccHnplex  processes  have  also  been  reported  |^9.  10,  llj. 

The  underlying  idea  of  d:/nariic  prograKsning,  first  introduced  1^  Bellman 
j  12  i  is  represented  hy  the  principle  of  optL-nality,  ^ich  states,  "An  optimal 
policy  has  the  property  that  whatever  the  initial  state  and  initial  decision 
are,  the  remaining  decisions  must  constitute  an  optimal  policy  with  regard  to 
the  state  resulting  from  the  first  decisicai"  1 12  |,  Qjr  means  of  this  principle 
and  the  method  of  interpolation,  a  multistage  decision  problem  is  replaced  by 
a  sequence  of  single  stage  decision  problems,  for  which  the  solution  is  drasti- 
cally simplified.  Using  the  methods  of  d^mamic  prc^ransdng,  Aris  has  solved 
a  variety  of  problems  concerned  with  optimal  cheraical  reactor  design  1  13j« 
Based  on  the  dynamic  prc^ramning  technique,  Ititten  and  Kemhauser  I  l^Jt  Aris 
j 15  I  and  Chien  j 16  j  have  presented  several  schemes  for  the  optLmization  of 
certain  types  of  topologically  complex  processes.  Aris,  Neahauser  and  'vllde 
1 17  I  presented  a  method  for  formulating  cyclic  and  brandling  processes  as 


initial  value,  final  value  or  two-point-boundary  value  problems.  Larson  j 18J 
has  developed  a  now-  approach  to  attack  the  probloms  of  continuous  character- 
istics, called  d:mamic  prograinning  with  continuous  independent  variables 
(state  IncroEient  ^jmanic  Prc^rainndng).  Ihis  procedure  has  a  high-speed  nffiiaory 
requirement  several  orders  of  magnitude  less  than  that  of  conventional  dynamic 
procrarining,  but  it  still  retains  the  general  applicability  and  other  desirable 
features  of  the  conventional  method. 

The  principal  purpose  of  this  thesis  is  to  show  how  the  raaxiimaa  pirinciple 
and  conventional  dynamic  programming  can  be  applied  to  topological3y  complex 
chffl^dcal  processes,  particularly  to  flow  type  chanical  reactor  systems.  T^ie 
study  is  divided  into  t^«>  parts.  The  first  part  is  an  investigation  of  optimal 
temperature  profiles  in  a  series  of  plug  flow  reactors  vri.th  either  fresh  feed 
to  or  intermediate  product  witlidrawal  from  an  internal  point.  A  first-order 
consecutive  irreversible  reacticsi  A  — *  B  ^  C  is  taking  place  in  the  re- 
actors. A  split  type,  two-point-boiaidary  value  problem,  ^rfiich  is  encoxmtered, 
is  solved  niimerically.  '0ie  second  part  deals  ^dth  the  combined  use  of  the 
raaximum  principle  and  dynamic  prc^rarardng.  General  working  equations  for  a 
sequential  union  of  these  two  techniques  are  developed.  'Kie  use  of  this  imion 
will  be  illustrated  by  ti-w  chemical  rate  processes  with  geometrically  complex 
features. 


2«   THE  MAXMJM  PRINGIPLB 

The  Baxiawia  prindixLo  was  first  proposed  ty  a  Russian  raathotaatician, 
Ponttyagin  111,  in  1956.  Pontiyagin's  original  method  «as  confined  to 
simile  cOTtinuous  processes,  that  is,  processes  without  any  geometrically 
ccxapleoc  structure  that  nay  be  described  by  a  systoji  of  first  order  differential 
equations.  Many  of  the  processes  encountoired  in  practice  are,  however, 
caaposed  of  several  inter-connected  branches.  Although  the  optimal  policy 
for  sudi  a  coRplax  process  can  soraetlwes  be  obtained  by  disecting  the 
process  into  several  simple  subprocesses,  independently  optimizing  each 
subprocesses,  and  thai  recorabining  these  into  the  overall  optimal  policy, 
it  is  more  desirable  to  use  an  algorithm  capable  of  obtaining  directly  the 
optimal  poli<^  for  the  entire  process  without  decomposing  the  process. 
The  following  essentially  follows  Fan  and  l-feng  i  8j, 

a»   Algorithm  for  Continuous  Complex  I^t-ocesses  I  8J 

Schanatic  representations  of  two  types  of  continuous  conpleEx  processes 
are  shown  in  Figures  1  and  2,  Figure  1  shows  a  ccmtinuous  conpleoc  process 
involving  a  separating  point,  that  is,  a  point  ^*lere  c«e  branch  of  the  path 
splits  into  several  brandies.  Figure  2  shows  a  complex  continuous  process 
involving  a  ccrabining  point.  A  combining  ix>int  is  a  point  t^ere  several 
branches  combine  into  one  single  branch.  Suppose  that  each  branch  in  a 
complex  process  is  described  by  a  performance  equation  of  the  following 
type 

dT'^i  <^—- V  h V  ^^ 

i  =  1,  2,..,,  s  ,   0  <  t  <  T 
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or  In  vector  foim 

^  =  f  (xje)  ,       0  <  t  <  T  (2) 

In  equation  (1),  s  and  r  are  the  dimensicais  of  the  state  vector  x  and  the 
decision  vector  0  respectively.  The  form  of  the  function,  f(x;0),  and  the 
langth  or  time,  may  differ  for  different  branches.  A  number  may  be  arbi- 
trarily assigned  to  each  branch  for  convenience.  Ihe  superscripts  In 
f^'^^Cx   I  e^^O  and  the  subscript  of  T.  represent  the  branch  number. 

The  arrows  at  the  junction  point  show  the  direction  of  the  flow  of  the 
process  streams.  Ihe  relationships  at  the  junction  points  beti/een  the  state 
vectors  of  different  branches  are  described  by  tlie  "jiinction  equations". 
These  are  (all  the  branch  numbers  are  referred  to  Figures  1  and  2)t 
(1)  Separating  points 

(3) 


(2)  Combining  points 

(n),^,    (n),  (1).  ,   (2)         (n-l)/n,   ^^      n.^ 
X   (0)  =  g   (x   (T^i^),  X   (T2) X    i\Jl))  W 

i^ere  the  function  g^"^(x^^^(Tj^))  is  the  vector  form  of  g("^(x|^^(Tj^), 

(k)         (k) 
xX     (Tj.)»«»»»  ^  fk)('^k^^'  ^^  "'^  ^  noted  that  the  dimensions  of  the  state 

s 
vector  may  be  different  fron  branch  to  branch.  Superscripts,  such  as  in 

s'^,  are  used  to  distinguish  t?»e  diffei*ences. 

In  a  complex  process,  there  may  be  sorae  branches  **iich  d<»i't  have 
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decision  acticai,  such  as  the  fresh  feed  stream  or  the  bypass  product  stream 
at  the  internal  point  of  a  series  of  tubular  reactors.  The  performance 
equation  for  such  a  non-decision  branch  is 

f  =  Hx)  (5) 

or  for  the  special  case  in  \i\i.ch  the  state  vector  remains  constant,  equation 
(5)  reduces  to 

i=o  «' 

A  ccsnplex  process  may  have  several  initial  points  and  final  points  and  the 
values  of  the  state  vectors  at  these  points  are  called  initial  states  and 
final  states  respectively. 

■Rie  optimiaation  problen  under  consideration  may  be  stated  as  follows: 
Given  all  of  the  performance  equations,  Action  equations,  lengths  (or 
tines),  and  initial  states  of  a  complex  process  c(»isisting  of  a  path  with 
N  branches,  b  initial  points,  and  f  final  points,  find  the  decision  vector 
function  of  each  branch  so  as  to  maximize  a  certain  linear  function  of  the 
final  states  of  the  process,  such  as 
f   s(^) 
te=l  i=l 

ytiere  the   summation  is  over  all  the  components  of  the  state  vectors  and  over 
all  the  final  states.  The  set  of  the  decision  vector  functiai  G^  (t), 
k  =  1,  2,».*,  N,  thus  found  is  called  the  optimal  policy.  Ihe  function  to 
be  maximized  is  called  the  objective  function. 

The  procedure  to  find  the  optimal  policy  for  a  coaplex  process  can  be 
divided  into  two  steps: 
Step  1.   Introduce  an  adjoint  vector,  z^  '(t),  and  a  Harailtonian  function 


t  7  4^^^h\). 


h'^'(x,  z,  e)  for  eac*  branch  satisfying t 

(k) 
^^^h..   e.  G)  =  I    4^k(^h4^^ -%l   0(^) 0%)       (7) 


1^1 


dz{^^(t)    aH^^^(x,z>G) 


.(k) 


IF— = ^;T^) •    i=l»2 s  (8) 

with  the  values  of  z^^   at  the  junction  points  and  final  points  satisfying 

(a)  separating  points 

(k)  .  (k) 
i   1   k=2  j=l   j^V-l;  3 

(b)  combining  points 

(n)  ^(n)  1=1,  2,,..,  s 

^i  ^^k^  "  ^,   rlET  ^J  ^°^  »       k  =  1.  2 (n-1) 

j=l  ^\ 

(c)  final  points 

1  =  12  s(^> 

«[^^(Tj^)  =  -  cf  >   ,  (11) 

k  =  If  2f*«*,  i 

Step  2.:     Find  the  decision  vector  function,  B^  \  from  the  following 
conditions 

H^     (x,2,e)  =  minlinura  at  eveiy  point  t,      k  =  1,  2,«.*,  N  (12) 

./If)  (Ir) 

In  other  words,  the  optimal  functions,  G^  ',  are  obtained  by  vaiying  9^  ' 
until  the  ralninun  value  of  H^^'  is  attained  at  each  point.  Ihe  proof  of 
this  algorithm  is  given  by  reference  18 J, 
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3,   APPLICATIONS  TO  FLOW  CHrStTCAL  H35ACTOR3 

It  is  known  that  for  a  sinjilo  reaction,  e.e.,  the  first  order  reversible 
reaction  A  < — *  B,  taking  place  in  a  plug  flo%i  reactor,  the  best  yield  is 
obtained  if  the  temperature  is  so  adjusted  as  to  raake  the  reaction  rate  as 
larce  as  possible  at  eacii  point  of  the  reactor.  This  is  the  so-called  disjoint 
property  of  tenporature  for  a  single  reaction  1 13  U     In  the  case  of  two  or  sjore 

competing  reactions  this  is  no  longer  true.  The  problom  of  a  first  order 

k     k 
irreversible  consecutive  reaction  \  — i^  "  ^C  has  been  extensively  dis- 
cussed in  the  %fork  of  Bilous  and  Ajmmdson  [  I9j»  They  derived  equations  to  be 
solved  simultaneously  vdth  reaction  onuati<»is  to  give  the  optimum  temperature 
profile.  They  use  two  methods,  the  first  based  on  an  axplicit  solution  of  the 
reaction  equation,  and  the  second  an  iraplicit  method  "which  may  be  applied  to 
the  reacticm  equations  thenselvos.  Both  involve,  however,  the  rather  difficult 
notions  of  functional  differentiation,  and  they  really  left  the  problem  un- 
solved. Aris  j 13  j  continued  the  work  by  the  nothod  of  dynamic  progranming. 
Using  the  dynamic  prograntning  technique,  he  obtained,  by  a  lengthy  derivation 
and  explanation,  a  set  of  characteristic  equations,  '<diich  can  be  easily 
obtained  by  direct  use  of  the  maxinum  principle.  Most  of  his  work  j  13j  was 
essentially  a  qualitative  analysis,  as  he  did  not  obtain  a  solution  of  the 
derived  equations.  Lee  i  20j  treat®3  this  pjroblen  by  the  conbined  us©  of  the 
maxinum  principle  and  a  naxiraura-seeking  nethod  cai  the  analog  coraputer.  In  his 
work,  the  rainLTaR  of  the  Hamiltonian  is  sought  directly  hy  a  random  search 
technique. 

All  of  the  ^wrk  dcscribod  above,  however,  ^-ras  confined  to  the  investi- 
gation of  the  optijJial  tomperatiu:'e  and  concentration  along  a  single  plug-flow 
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reactor.  If  there  are  two  plug  flow  reactors  arranged  In  sertes  with 
product  withdrawal  from  the  first  reactor,  or  firesh  feed  into  the  second 
reactor,  the  optimization  problem  canH  be  treated  routinely*  This 
geometrically  complex  arrangfflnwit  of  reactors  will  render  the  problem  more 
complicated  and  difficiilt,  especially  with  respect  to  the  numerical  solution. 
The  reactor  arrangeraents  (or  models)  encountered  in  this  investigation  may 
often  be  found  in  industry.  It  is  often  desirable,  for  econanic  or  practical 
reasons,  to  use  two  reactors  in  series,  each  with  a  smaller  reactor  volume 
instead  of  using  one  large  reactor.  Also,  vhen  ve  need  a  certain  amount 
of  product  with  lower  conversion,  or  ^en  we  have  two  sources  of  raw  mate- 
rials, different  both  in  ccmposition  and  availability,  we  might  have  to  set 
up  the  reactors  in  the  way  shown  in  either  Figure  3  or  arranging  two  tubiilar 
reactors  in  series  ■^dth  fresh  feed  fed  to  an  internal  point. 

a»   Two  Tabular  Flow  Reactors  in  Series  with  IDitenaediate  Product  withdrawal. 

C<aisider  a  complex  continuous  process  consisting  of  two  tubular  flow 
chemical  reactors  connected  in  series  with  exit  flow  frcan  an  intermediate 
point,  as  shown  in  Figure  3»  The  reactions  taking  place  in  the  reactors 

are  considered  to  be  first  order  irreversible  consecutive  reactions  of  the 

k,      kp 
type  A *■  B     »  C.  Ihe  length  of  reactors,  L^  and  Lg,  and  the  cross- 

sectlcwal  areas  A^  '  and  A^  '  are  all  specified,  that  is,  the  volumes  of 
the  two  tubular  reactors  are  fixed,  or  alternatively  the  holding  time  of 
each  reactor  is  fixed  for  specified  volumetric  flow  rates  of  feed  and  inter- 
mediate product  streams.  It  is  desired  to  maximize  the  total  amount  of 
product  B  from  the  intermediate  and  end  product  streams.  An  equivalent 
objective  functicm  is  to  maximize  the  concentration  of  B  at  the  ends  of 
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both  reactors*  since  the  volumetric  flow  rate  of  each  stream  is  fixed,  Ihe 
flows  In  the  reactors  are  assumed  to  be  plug  flows,  i.e.,  it  is  assumed 
that  a  velocity  gradlaat  in  the  radial  direction  and  diffusion  in  the  radial 
and  in  the  axial  directions  do  not  exist  in  the  reactors. 

■Rie  steady  state  differential  material  balance  for  reactant  i  in  the 
differential  section  de  of  the  kth  branch,  k  =  1,  2,  can  be  written  as 
Rate  of  input  of  component  i  to  the  differential  section  dj, 
W  the  kth  branch  "' 

Rate  of  output  of  comp<»jent  i  frcai  the  differential  section 
"  W  of  the  kth  branch  •' 

Rate  of  productioi  of  ccmponent  i  in  the  differential  section 
^di  of  the  kth  branch  ^ 

=  0 
Ihe  derived  equation  is 


or 


-37—  =  ^  Rf^   .     i  «  1,  2,  and  k  =  1.  2  (12) 


xi^ere 


Oc) 

x^  '  =»  molar  concentration  of  reactant  i  in  the  kth  branch 

fk) 

V*   a  volumetric  flow  rate  in  the  kth  branch 

(k) 
A^  '  =  cross  sectional  area  of  plug  flow  reactor  in  the  kth  branch 

(k) 
R^  '  =  rate  of  producticai  of  caapcment  i  in  the  unit  volume 


v^ 


I  a  axial  position  in  the  plug  flow  reactor  (measured  fran  the  inlet) 
If  X,  and  Xg  are  the  molar  concentrations  of  A  aiKi  B  respectively, 
the  perforraance  equaticais  for  eacdi  branch  can  be  written,  according  to 
equation  (12),  as 

"dr^  =  n:iTt^  J  ^^^ 


.(1) 


^.a":^.JC1)l  (14) 

^r^  -  TOT  1.^  J 


-IT- =  7^1%  ;  ^^^^ 


(2) 


^2     a(2)^,2) 


Since  the  mechanism  of  reaction  is  the  same  for  both  branches  1  and  2,  on© 
can  write 

4«  =  -  fc?=>  ^'>  (17) 

gCk)  .  i^k)  ^k)  .  ^k)  j^k)  _     k  =  1.  2  (18) 

\#iere  k^  and  kg  are  tenperature  dep«ident  specific  reaction  rate  constants 
of  the  Arz^ienious  type  ^iven  as 

\  % 

k^  =  kj^Q  e    ^^  ,  ^2  =  ^20®^  (19) 


lAere 


k.Q  =!  frequency  factor*  hr 

E.     a»  activaticxi  energy,  cal/raole 

6   =  temperature,  *^ 

R   =  gas  constant,  cal/mole  ^K 
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If  the  holding  time  ts  defined  as  t  =  ^,  substitution  of  equations 

(17)  and  (18)  into  equations  (13)  through  (16)  yields 

(20) 

^^  _  (1)  (1)  ^  Al)     (1)     JD^QS   ^ 

(21) 

^  (2) 

(22) 

^=.f:4^'-^^)4^> 

(23) 

where  o-,  and  o-g  are  the  inlet  concentration  of  A  and  B  respectively. 

Since  no  transf ornatian  is  taking  place  and  no  control  action  is 

imposed  on  branch  3,  the  perforaance  eqixations  for  that  branch  are 

dxp>   dx^3) 
di  -  dt   "° 

(24) 

Since  the  outlet  stream  of  cos^nent  i  from  reactor  1  splits  into  two 

strearas  at  the  separating  point,  the  concentration  of  compon«it  i  in  both 

streans  at  this  junction  must  be  equal,  that  is, 

*[^^T,)  =  xf  >(0) 

(25) 

x[^^(T3^)«x[^^(0) 

(26) 

4^\t{)  =  4^\o) 

(27) 

4^\h)  =  4^\o) 

(28) 

lAere 

T   -     "^ 

^  "  v<^)/  a(^) 
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is  the  total  holding  time  of  the  kth  reactor,  k  ==  1,  2. 

According  to  ©quatiwis  (7)  and  (8),  the  HawiltOTian  functicai  and  the 

adjoint  variable  functions  can  be  xTrltten  as 

„(1)    .  (1)  (1)  (1)  ,  ..  (1)  (1)   .  (1)  (1),  (1) 

(29) 

H<^>  =  -  .f  xf  .f  +  (kf  xf  >  -  k<^'x|^')  4'^ 

(30) 

H«^0 

(31) 

#=^^M^)-.p>4" 

(32) 

dt   =  "2  "2 

(33) 

(2) 

(3^) 

(35) 

dt    "   dt    =° 

(36) 

Accomiing  to  equation  (9),  the  adjoint  variables  at  the  separating  point  are 

givwi  by 

(37) 

Since  it  is  desired  to  maximize  the  intermediate  product  B  at  the  ends  of 

branches  2  and  3»  the  objective  functicm  is 

5  ^  4^\'^z) + 4^\t^)        I.  ■/'-'■'" 

■  *'     ■    y  ■  .i  '^ 
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Consequently 

of  >  =  =p>  =  1 

(38) 

42)  ,  c<3>  .  0 

(39) 

Substitutioa  of  equaticais  (38)  and  (39)  into  equatloi  (11)  yields 

at  the  final  points 

4^^(T2)  =  4^\t^)  =  0 

(^0) 

8^2^  (Tg)  =  4^^(13)  =  -  1 

(^1) 

The  optimal  ten^rature,  0(t),  is  determined  fraa  the  conditicm  that  H 

is  lainicaum  at  every  point  of  each  branch.     Assuming  that  the  miniimuii  of  H 

occurs  at  the  staticmary  point,  the  optiuial  decision  for  branch  1  can 

be 

obtained  from 

m^^^       ^a),(i)%     ,  .%      (1)     ^     ,(!)> 

4"=- 

0 

or 

<4'M'^  -  .f  M'^)  ^0^1  -p  (-v^^'^) 

-  z^^^xf  ^  k2o':2  ^  (-Sg/RQ^^b  »  0 

(^2) 

Solution  of  this  equation  for  0^  '  gives 

^(l)/xx         A""^x                          1 

^2     ^2 

(^3) 

-uhere 

*^20^'2 

Id 


Slrallarly,  for  branch  2,  one  obtains 


Substitution  of  equation  (^3)  Into  equations  (20),   (21),   (32),  and  (33), 
and  eqiiaticai  (^)  into  equations  (22),   (23),   (34),  and  (35)  yields 

-dt~  =  -  ^10  L^  "07  ^ — aT~^J    ^1  ^"^^^ 

^  ^^2 

d,(l)  ^(1)    Jl)  .  Jl)  x_ 

^  U        Ir^         ^"^2        -'l      .h      (1) 

(1)     g(l)  .     (1)    X 


dji)  3J1)    2(1)  ,  5.(1)    X 

-i~  =  '^o  L?  737  ^  '   ,(1)     )J  '  ^4'^  -  4'^>  C^7) 


.  (1)  ^(1)     (1)      (1)    . 

dJ2)  (2)     ^(2)  (2)     j^ 

^  =  -  ^0  L^  227  ('     ,(2)       )J  '  4'^  (^9) 

^  *2 

^(2)  ,(2)        (2)  (2)     . 

^      _  ^      ,  -  *1       /f2_-2L.^  1^1  ,(2) 

*2  ^2 

(2)     ^(2)  .     (2)     j^ 

-  ^0  L?  -W  <"       (2)        )j  '  4'^  C50) 

*2  *2 
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(2)                      (2)       (2) 

OZ-f                                       Xy             Zp 

6.i       ho  'J  J2)  ^ 

-j#,/i  (4^)  -  4^h        (51) 

^2 

.   (2)                       (2)       (2) 

"2 

viiere 

,      'i 

X    =     '^     . 

^l-E^.E, 

2      '^-'^ 

For  this  ostein  of  eigjit  difforoitlal  eqtiations,  el^t  boundary  conditiais 

are  needed,  ^^idi  can  be  suunnarizod 

as  follows: 

x^^ho)^a^ 

(20) 

4^\0)=a^ 

(21) 

x[^^(T^)  =  x[2)(0) 

(25) 

7:^^h\)  ^  4^ho) 

(26) 

^^^(Tg)  =  0 

(40) 

zf  )(T2)  =  -  1 

(41) 

zf^hl^)  =  z^^^CO) 

(53) 

z^d;^)  =  z|2^(0)  -1 

(54) 

Hquations  (53)  and  (5^)  are  obtained  by  conbining  equations  06),   (40),  and 

(41)  vith  equation  (37). 

The  nuraertcal  calculations  for  the  system  described  by  equaticais  (45) 

through  (52)  are  illustrated  idth  the  use  of  the  following  datai 

20 


^10  =  °'^55  *  ^^^   (minute)'^ 
kgQ  =  0.4^  X  10   (rainute) 
E,  =  18,000  cal/mole 
EL  »  30,000  cal/mole 

R   =2  caX/mdXe   °K 

T,  =  10  minute 
1 

!«  =  8  minute 

a_  =0.72  moles/liter 

tfj,  ■  0#2t*  moles/liter 

An  ai^roxiraate  iterative  Runge-Kutta  method  was  used  together  ^dth  a 
backward  integration  procedure  to  obtain  the  numerical  results.  The 
nvuaerical  procedure  is  as  follows: 

1.  Differentiate  equation  {k5)   with  respect  to  x£  ^(T^^),  equation 
(46)  with  respect  to  x^Ht^),   equation  (49)  with  respect  to  x£^'(T2)f  and 
equatiwi  (50)  with  respect  to  x^^  (Tg)  to  obtain  four  additional  first 
order  differential  equations.  They  are: 

.^(1)  .M^\  Jl)    Jl)   .  ,(1)  X 

"^^  °  THTTT  =  "  ^0  ^^1+ ^H^  75T  ^ lir-^J  ^1 

dxi  (Ti)  Xg      Z2 


(5^a) 


._(1)   w^  .  fl)   (1)    (1)  >„(1) 

dXj  ^T^;  Xg      Z2      ^ 


^(1)   (1)    (1)  . 
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dC 


(2) 


dxf> 


(2)        (2)       _(2) 


^2) 


(2)        (2)  (2)     ^.„(2) 


dt 


dt 


dx^  ^(Tg) 


"  no  ^1  L^  I2T  ^ 


r<^ 


12) ^J  ^727^  °2 

*2  ^ 


(2) 


(5^c) 


(2)     ^(2)  ^  ^(2)     j^ 

+  ^20   (^2  -  1>  L^  ^  ('       (2)^     )J  '  4'^ 

Xg  «2 


(5^) 


%*iere 


.(2)  _  J^ 


(2) 


dx^^^ 
(1)  2 

°2    =  — m — 

2     d^^Pcy 


G 


dx^"'(T2) 
The  boundary  conditions  for  g|-^',  G^',  G^ ',  and  G^)  are 


XD 


t=T, 


dx. 


(2) 


T^^v  =  r(^ 


dx^^^'T2) 


t=T« 


dx^2) 

»  1.  G^  ^(Tg)  =  -75^ 

2    2    dxf^CTg) 


(5^e) 


=  1 


t=T, 


(5^f) 


t=T. 


Thus,  for  these  four  additional  equationsi  equations  (5^a)  through  i^^)t  one 
has  four  additional  boundary  conditions,  equaticai  (5+f)»  In  the  derivation 
of  these  equations,  equaticms  (5^a)  through  (5^),  the  assumption  that  oily 
x£^'  is  a  function  of  x^^^(T^),  only  x^"'''  is  a  function  of  xP^d^),  only 
x^2/  is  a  function  of  x.[  (Tg).  and  only  t^'   is  a  function  of  xi^^CTg)  is  made. 
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2,      Assume  values  for  x^^^CTg)  and  x^^^CTg).    Note  that 


and 


4^^(T2)  =  0  . 


o[^\tz>  =^  • 


.(2). 


z^-^CT^)  =  .  1 


G^^^CTg) 


3,      Use  the  foixrth-order  Rimge-Kutta  method  integrating  backward 

(2) 
equations  (^9)  throu^  (52),  equations  (5hj)  and  (5W[)  to  obtain  Xj^  (0), 

x|2)(0).  4?>(0),  z|2)(0).  g(2)(0).  and  g(2)(0). 

if.   Obtain  x^-^^Tj^),  4'''^^^1^'  ^''''^^l^'  ^^  4^^^%^  ^^'^  equations 
(25),  (26),  (52).  and  (53)  respectively.  Note  that 

.(1), 


G^^^(T^) 


and 


'2  (V 


5,  Ihtegrate  backward  anplqying  the  fourth-order  Runge-Kutta  method 
to  obtain  x^^>(0).  x^iO),  4^>(0).  4^>(0).  C-J[1^(0),  and  g|^>(0). 

6,  Use  the  derivative  aRoroxitaations 


0<"(0) 


«» 


g(i)/o)--4S— 


x^^^(O)  -c^i 

Jl)..  ^       ,(1) 
^     (^l)  -  ?i 

t  =  0        ^2     ''^l^       -2 


fa 


to  compute  gP  and  ?i^^  liiich  are  better  approxijnations  to  the  values  of 
x}^\t,)  and  x^  (\)  ^^^  those  obtained  in  Step  4. 

7.  (btain  froi  equations  (25)  and  (2?)  closer  approximations  to 

the  correct  values  of  x£^^(0)  and  x^^  (0).  Lot  these  values  be  f£  ''  and  fg 
(These  values  are  identical  with  those  obtained  in  Step  6). 

8,  Use  Uie  modified  derivative  a^-proximatlons 


D 


44  (Tg) 


0    *2  ^^2'   ^2 


to  dbt&in  5P'  and  5^^^  whloh  are  the  new  trial  values  for  x^^'dg)  and 

9»   Repeat  steps  3  throu^  8  by  using  the  new  trial  values  xmtil 
the  c<sigxited  values  of  xJ^(O)  and  xi'^'(O)  are  appraximately  equa-T  to  er- 
and  a,  respectively. 

10.  Calculate  the  optiraal  decisions  9   (t)  and  8   (t)  according 
to  equations  (43)  and  (^). 

This  oonputational  jwooedure  was  prograjmaed  for  an  IBM  1620.  "Cie  main 

diffictilty  of  the  ccsaputation  is  the  guess  of  the  first  assumed  values  for 

(2)        (2) 
3^    (^2)  and  X2  (Tg)*  Because  Meaningless  results  \jould  be  obtained  if 

either  x^t  Xg*  ^"^  ^^''^^^^z   -^^'''*  negative  in  any  branch,  and  a  very  long 
time  would  be  needed  to  obtain  convergence  to  the  desired  values  if  tae 
guess  was  too  far  away  flrcan  the  correct  values,  a  careful  choice  of  the 
value  of  X'j  '(T2)  and  x^  '(T2)  was  necessary.  To  do  that,  a  table  of  optimal 
solutions  without  any  iteration  yoiS  obtained  first.  Ihis  was  done  as  follows t 

Artitrarily  but  reasonably  pick  sets  of  xi  '(Tg)  and  xi  '(Tg)*  For 
each  set  of  x^^'CTg)  and  xl^^CTg)  together  with  given  values  of  s^^^CTg) 
and  8^  '  (Tg)  employ  the  fourth-order  Runge-Kutta  nethod  and  backward  inte- 
gration to  int^rate  equations  (^5)  through  (52)  to  obtain  a  set  of  values 
of  Xj^   '(0)  and  Xg  (0).  Each  backward  integration  gives  the  optimal  path 
correspwiding  to  the  initial  conditicms  thus  computed. 


2f¥ 


.(2)/m_\  „J  „(2), 


Iht^iratinc  in  this  way  for  many  sets  of  x£  '(Tg)  and  x^  (Tg)*  an 
optimal  soluticai  table  is  obtained.  It  is  worthviiile  to  raonticai  hei^  that 
a  slight  change  of  either  x^^CTg)  or  xi^CTg)  would  result  in  a  large 
change  in  x|^(0)  and  x|  '(0). 

Ejy  investigating  the  above  mentioned  table  and  noting  the  initial 

conditions,  x|-^^(0)  and  x|^^(0),  the  values  of  x^Ht^)   and  x^^^CTg)  for 

(Z)  (2) 

the  iteration  calculation  were  chosen.  -^Veani  though  x^^     (Tp)  and  xi  (Tp) 

were  selected  in  this  manner,  the  convergence  of  the  calculation  still  was 

such  that  the  answer  could  not  be  obtained  in  a  reasonable  period  of  time. 

To  overccBJie  this  difficulty,  a  modified  derivative  approximation  for  the 

(2)  (2) 

values  of  Gj^  '(0)  and  Gg  (0),  i^ich  is  given  in  Step  8  of  the  numerical 

procedures,  was  used.  The  reason  that  this  modification  t«s  needed  mi^t 
be  due  to  the  fact  that  there  are  jimction  equations  relating  the  end  values 
of  the  state  variables  and  adjoint  variables  of  branch  1  xdth  the  initial 

values  of  branch  two.  In  each  trial  calculation,  new  trial  values  of 

(2)        (2) 
*L  ^^^2^  ^^"^  ^  ^^2^  ^^^^  obtained  by  the  equations  in  Step  8,  i^iich  involve 

the  state  variables  of  both  branches.  It  could  be  because  of  this  "sudden 

Ju^p"  or  "smoothlessness"  v^ich  made  the  convergence  difficult  to  obtain. 

B^  employing  the  modified  derivative  approximations,  about  six  minutes  of 

computer  time  \ns  consumed  for  each  trial  calculation.  Ihe  rate  of  ccmver* 

genoe  became  slower  as  the  computed  values  of  x|  '(0)  arei  xi  '(0)  converged 

toward  the  givoi  initial  conditions.  During  the  calculation  new  trial 

values  can  also  be  put  in  to  make  convergence  faster  by  observing  the  trends 

of  the  results  coming  out.  Ihe  results  of  computation  are  presented  in 

diraensionless  time,  that  is 

VTj. 
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as  sho«Q  in  Table  1  and  Figures  >*  and  5«  For  the  first  assumed  values, 

JcJ^^CT,)  «  0»115,  x^^^(T  )  =  0,72,  the  coapater  took  about  tuo  hours  to 
X        £■  2    2 

obtain  the  solution  shown  in  Table  1« 

Figure  5  shows  that  the  optlnal  teog)erature  for  reaotor  2  is  almost 
constant  along  the  reactor.  For  design  purposes,  this  indicates  that  iso- 
thermal operaticai  in  the  second  reactor  will  almost  raeet  the  optimum  condi- 
tions* Even  though  this  fact  depends  on  the  specific  nuraerical  data  chosw 
and  might  not  always  be  true,  one  still  can  safely  assume  that  the  optimal 
temperature  profile  in  the  sec<aid  reactor  for  this  kind  of  systan  will  pr<^ably 
always  be  flatter  than  that  of  the  first  one. 

rt  should  also  be  noted  that  the  geometry  of  the  reactor  arrangement  in 
this  gystea  is  oquiTalent  to  (xae  plug  How  reactor  of  total  holding  tine 
T  »  \f  ^2*   with  a  side  exit  stream  at  t  =  T^.  The  aystoa  has  the  character- 
istics of  a  plug-flow  reactor,  providod  that  the  product  c«icentratic«i  fran 
the  exit  stream  is  not  to  be  maximized*  If,  as  in  this  problem,  the  linear 
eoaift)inatlon  of  end  product  concentrations  from  branches  2  and  3  is  to  be 
naxiraiaed,  such  an  equivalence  no  longer  exists.  This  fact  can  be  visualized 
froia  the  boundary  conditions  given  in  equations  (^1)  and  (^).  If  only  the 
Old  product  concentration  from  branch  2  is  to  be  maximized,  the  condition  in 
equation  (41)  should  be  changed  to  zp'(T-)  =  0,  \Aioh  coBibines  with  equation 
(37)  to  give  Bg  ^(Ti)  ^-  b|^'(0)  rather  than  the  value  given  in  equation  (,5k). 
This  changes  the  problem  frcm  that  of  a  simple  exit,  because  both  the  optimal 
concentration  and  toaperature  are  affected  by  the  values  of  the  adjoint 
variables. 

b«   TWO  Tubular  Flow  Reactors  in  Series  with  Intermediate  Fresh  Feed 


a( 

Table  1. 

Computational  results  for  the  reaction  A-»-B-»C  taking  place  in 

tw5  tubular  reactors  in  series  with  Interrwdiate  product  vdthdrawal. 

4^) 

4==) 

4« 

f) 

9(2) 

1.0 

0.11315 

0,70630 

0,00000 

-1.00000 

327.66 

0,9 

0.11901 

0.70387 

-0.04917 

-0.99515 

327.72 

0*8 

0.12522 

0.70111 

•0.09591 

-0.99027 

327.79 

0.7 

0.13178 

0.69801 

-0.14033 

-0.98538 

327.35 

0.6 

0.13873 

0.69455 

-0.18253 

-0.98046 

327.93 

0.5 

0.1i^609 

0.69069 

-0.22^61 

-0.97552 

328.00 

O^k 

0,15389 

0.68640 

-O.26O67 

-0.97055 

328.08 

0.3 

0.16217 

0,68166 

-0.29679 

-0.96555 

328.17 

0.2 

0.17097 

0.67642 

-0.33105 

-0.96051 

328.26 

0.1 

0,18031 

0,67064 

-0.36355 

-0.955^^2 

328.36 

0.0 

0.19026 

0,66429 

-0.39^36 

-0,95030 

328.47 

JX) 
*2 

J^^ 
^ 

*2 

e<i> 

1.0 

0.19026 

0.66429 

-0.39436 

-1.95030 

33^.15 

0.9 

o.;?i205 

0.65235 

-0.55271 

-1.92134 

334.39 

0.8 

0.23688 

0.63752 

-0.69462 

-1,89178 

334.68 

0.7 

0.26538 

0.61919 

-0.82149 

-1,86144 

335.02 

0.6 

0.29831 

0.59658 

-0.93455 

-1,83009 

335.42 

0.5 

0.33678 

0.56863 

-1,03493 

-1.79738 

335.92 

0.^ 

0.38229 

0.53383 

-1,12361 

-1.76282 

336.54 

0.3 

0.43717 

0.48988 

-1,20146 

-1.72561 

337.35 

0.? 

0.50520 

0.43303 

-1,26921 

-1.68435 

338.48 

0.1 

0.59363 

0.35607 

-1.32739 

-I.636I6 

340.20 

0.0 

0.72084 

0,24071 

-1.37616 

-1.57259 

343.58 

» 
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In  Sectloi  a*  a  continuous  eon^plex  process  involving  a  separating  point 
has  been  discussed.  In  this  section,  a  similar  problem  will  be  considered, 
i.e.  a  reaction  A  -—►  3  —^  G  is  taking  place  in  the  reactors,  but  the 
geometrical  arrangement  involves  a  ocxabining  point  rather  than  a  separating 
point.  A  schematic  diagram  for  this  process  is  shown  in  Figure  6.  In  this 
figure,  v^^'  represents  the  volumetric  floxir  rate  through  the  kth  branch.  It 
la  seen  that  the  reactor  arrangeraGnt  is  similar  to  that  of  Socticai  a,  except 
an  additional  fresh  feed  streams  is  introduced  at  an  intermediate  point  of 
the  two  reactors  rather  than  a  ^-dthdrawal  of  side  stream  pi*oduct.  The 
lengths  of  the  reactors,  L,  and  L,  and  the  cross-sectional  areas  A^  '  and 
A  ■''  are  all  preassigned,  that  is,  the  volumes  of  the  reactors  are  fixed.  In 
this  study,  v^-'*^  and  v^-^',  tlie  volumetric  through-put  to  the  first  reactor 
(branch  1)  and  the  second  reactor  (branch  3)  respectively,  are  fixed,  i.e., 
the  total  holding  time  in  both  reactors  has  been  specified.  It  is  desired 
to  maximize  the  intermediate  product  3  at  the  end  of  the  second  reactor  by 
proper  choice  of  toaperature  along  the  ti-ao  two  tubular  reactors  \jith.  given 
inlet  concentrations  of  both  feed  sti^ams.  The  floira  in  the  reactors  are 
again  assumed  to  be  of  the  pliig  flow  type* 

From  the  material  balance  of  each  component,  as  carried  out  in  Sectl<»i  a, 
the  performance  equations  for  branches  1  and  3  liiere  the  reactions  take  place 
can  be  T-n:^tten  as 

^2^  =  -  k^l^x^^^  (55) 

dt     ^1  '^l  ^^^' 

^2  .  (1)  (1)   .  (1)  (1)  /-,x 

-dT-  ==  ^1  '^l  -  4  4  (56) 

H,(3) 
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^-  =  J3)x,0)  .  kp)  J3)  (58) 


vith  the  given  initial  concentration  of  A  and  D  as  x£  ^(0)  =  or^^,  x^  '(0)  =  org 
respectively!  x,,  x^,  t,   k^  axvd  kg  are  the  sane  as  defined  in  Section  a. 
Since  there  are  no  transformations  taking  place  and  no  control  actions  imposed 
in  branch  2,  the  porfomance  equation  for  branch  2  is 

^r^    dt        • 

t^ere  P^^  and  Pg  ^tre  the  concent  rat  ic«s  of  A  and  B  respectively  in  the  inter- 
mediate feed  and  are  generally  different  fron  the  values  of  o-^  and  c^g*  The 
juncticai  equations  at  the  coiabining  point  are  derived  based  on  the  principle 
of  ccxiservation  of  nass.  The  application  of  this  principle  (the  material 
balance)  at  the  combining  point  yields 

v«3)43)(o)  „  ,(«,(!) (,^)  +  v(2)x(2)(T2)  (60) 

liiere  v^^s  k  =  1,  2,  3  is  the  volumetric  flow  rate  through  the  kth  branch. 
If  we  let  Yi  »  v^  Vv^^'  and  7,  =•  v  /v'^',  equatiwi  (60)  can  be  rewritten 
as 

x^^^  (0)  «  Yix}^^  (Ti)  +  Y^xf  ^  (Tg)  (61) 

According  to  ©riuation  (59)  the  concwitraticsis  of  reactants  1  and  2  (A  and  H) 
remain  constant  along  branch  2,  It  can,  therefore,  be  replaced  by  a  acre 
general  forra  as 

x[^^0)  =  xi^ht)  «  xf  ^(Tg)  =  3;|,  (62) 

Substituting  equation  (62)  into  equation  (61)  yields  a  junction  equation 
at  the  combining  point  for  reactant  1  as 
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4^\o)   =  Yl==l"(Ti)+  VjSl 

(63) 

Slrailarly,  for  roactant  2  one  obtains 

4^^(0)  =  Y^xf^(Ti)4-V232 

(6^) 

The  Initial  states  are 

xj[^^0)»tti 

4^ho)»a^ 

(65) 

Accordinc  to  equations  (7)  and  (8),  the  Hajniltonian  functions  and  the  adjoint 

variable  functions  can  be  i-nritten  as 

H<^>  =  4^>(-lcf^^))+.f) 

(kf^: 

=^1 

-4" 

4^^) 

(66) 

H<^)  =  0 

(67) 

H<3)„43)(.,(3),p))^^(3) 

(kp^: 

'=1 

-.p> 

4'^) 

(68) 

di   "  n  '^l    "^2  ^2 

(69) 

dt    ^2  ^2 

(70) 

(71) 

^  =  k53U3)  .  k^3)  (3) 
(ii   ^  ^   *^2  ^2 

(72) 

^-?'4^> 

(73) 

According  to  equations  (10),  (63).  and  (6^^), 

one 

has  at  the  coriining  point 
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4^^(Ti)=Yi4^^0)  (75) 

sine©  the  problem  is  to  raaximlze  the  concentration  of  product  B  at  the 
outlet  of  branch  3,  the  objective  functicai  is  x^-'^iTj),     Consequently 
ci^'  =  0  and  cP '  =  1.  Thus,  according  to  equation  (11) t  one  has  at  the 
final  point 

4^^T3)  =  0  (76) 

4^^(13)  »  -  1  (77) 

The  optimal  temperature,  8(t),  is  deterrained  in  the  saiue  manner  as  is  givoi 
in  Section  a*  The  result  is 


^""^^'  ■  '^^'      4^)(4^)  -  4^)) 


^L^    50j(i0  ^J 

^ere 

^0^1 
k  =  1,  3   and   5  =  r — rr 

*^20''2 
If  this  result  is  substituted  into  equations  (55)   throu^  (58)  and 
equations  (69)  thirou^  (73),  one  obtains 


rfjl)         J^^    ,(^^   ,^^^  X   .  ^ 

-dr-=-^oL^;757^ — 7ir""^^J   "^ 
.ji)  ^(1)  ^(1) .  -d)  X    , , 

^    ^   ,  ,  ^   r!2Jlfl_N  1^1  ,(1) 


(79) 


(1)  ^(1)  .  (1)  3^ 

-  ^20  L^  TiT  ^'  (1)  ^J  '  4'^   (3°) 


dx.(3)        ~(3)  „(3)  _  J3) 


il-..  , , ^1  ("2  : '1  ) h „(3)  (31) 


-dT- =  ^20  L^  TIT  ^ — Z^r^^i   =^2  ^^^ 

.J3)  J3)   (3)   (3)  . 

^       *2 

-^  =  ^20  L^  %T  ^-^-75r-^J   4'^  <^^) 

^       ^2 

To  solve  this  system  of  eight  differential  equations,  one  needs  eight 
bo^jndary  conditions  x^icJ-i  are  given  by  equations  (63) »  (6if),  and  (65)   and 
equatlOTis  (7^)  through  (77)  and  are  aunnarized  as  foHovras 

4-^^(0)  =  Yixf^(Ti)+.  Y23l  <^3) 

4^^(0)  =  vi4^^(Ti)+  Y2&2  ^^^ 


4^^(0)=c.^ 
4^^(0)=c.2 


(65) 
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^V=Yl4^^0)                             (74) 

4^^ 

>(T^)  =  Yi4"^^(0)                            (75) 

^^■ 

\l^)  -  0                                   (76) 

h-r^)  =  -  1                             (77) 

The  nnmertcal  calculations  for  the  system  described  by  equaticMs  (78) 

through  (86)  i 

ire  illustrated  by  using  the  following  datat 

ho 

«  0,535  X  10^  minute"^ 

^20 

»  0,^+61  X  10-'^  minute"^ 

h 

=  13,000      cal/nole 

h 

=  30,000      oal/nole 

R 

-  2          cal/mole  **K 

h 

=  8          minute 

h 

=  10         minute 

"1 

=  0.72        moles/liter 

"z 

a  0,2i^        moles/liter 

h 

=  0,7        moles/liter 

h 

=  0,25       moles/liter 

An  approximate  iterative  Runge-Kutta  method  together  with  a  backward 

integration  procedure  \m.s  used  to  obtain  the  nvuaerical  solutions.  The 

numerical  procedures  are  similar  to  those  used  in  Section  a,  except  that 

in  obtaining  : 

t[-^^(T^).  x^-^-^dj^),  z^Ht^),   and  z^h\)   one  has  to  use  the 

boundary  conditions,  equations  (63)»  i6^)»   (7^)  and  (75)   respectively  (see 
Stop  k  of  Section  a).  These  boundary  ccmditions  render  the  computation 
nuch  more  difficult.  As  one  can  see,  equations  (63)  and  (6'f)  give 

x^^Ht^)  =  i.  (xJ^^(O)  -  Yg^i)   .  ^37) 

x^l>(T,)  =  i  (4^>(0)  -  Y2P2)   .  ^^^ 

^1 
rather  than 

4^hTi)  ^  x[^\o)  (25) 

4^^(11)  =  4^^  (0)  (27) 

A  slight  change  in  x^\t:J)   and  xp^(T-)  would  greatly  affect  x^-'-^(Tj^)  and 
Xg  (T,)  tiiich  irere  responsible  for  the  later  calculations  for  branch  1. 
I'feaningless  results  might  result  due  to  the  appearance  of  negative  values 
of  X,,  Xpi  or  (zp  -  24)/22»  Ihe  most  critical  part  of  the  calculation  was 
uh«n  the  corapiter  proceeded  fr(»x  the  calculation  of  branch  3  to  branch  1» 
A  raore  elaborate  table  of  optinal  solutions  as  moitioned  in  Section  a  ivas 
needed  to  make  a  suitable  choice  for  the  first  assumed  values  of  xV'(T^)  and 
Xg  (T«)  for  the  final  iteration  calculaticm,  Which  employed  the  approximate 
iterative  Runge-Kutta  method. 

In  obtaining  the  new  trial  values  for  xj^^   (T^)  and  x^^d^)  for  the 
iterati<»i  calculation,  it  was  necessaiy  to  modify  Steps  6,  7i  and  8  as 
follows:  (note  that  one  has  branch  3  instead  of  branch  2  for  representing 
the  second  reactor) 

6'   Use  the  derivative  approximatiwi 


of  ^(0)  =  -^S 

dxf^(Ti) 


xf^O).a^ 

"*  J1)(T  ).,(!) 
t  =  0   ^  ^^1^  h 
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oa)(o)  =  -^ 


to  compute  ^^^   and  F^-*-^  v;hich  are  better  approxiraticjns  of  the  values  of 


x^^^(T;^)  and4^^(T;^) 

7»   Obtain  from  equations  (63)  and  (6^)  closer  approxiraations  to 
the  correct  values  of  xp^(O)  and  x^^'(O),  i.e., 


^'^y^^Ky^i 


.p> 


Yi§f  ^  Y?32 


(89) 
(90) 


8*   Use  the  modified  derivative  approxiraations 


0,(3)  (0)  _  ^ 


Gp^O)  = 


.(3) 


1  ^^3^ 

dxp) 
dx^^Tj) 
.(3) 


|(x^^^(0)  -TI^^O 
t  =  0    ^1  ^^3''   -1 

|(xp^(0)  -  T]pb 
t  =  0    ^  ^^3^   '2 


(3). 


to  obtain  ?|-^'  and  fo  »  ^^°^  *'*®  "^^  "®w  trial  values  for  x^-^'d^)  and 

As  can  be  seai,  duiring  the  calculation  there  xms  another  "sudden- jtu^" 
in  Steps  7'  and  8«.  "Riis  nd^t  also  be  one  of  the  reasons  idiy  the  computation 
was  so  unstable  and  so  slow  to  convorge* 

This  eoD^tation  was  also  prograitimai  on  an  I3M  1620  conputer.  About 
sevon  minutes  were  consumed  for  each  trial  calculaticsi.  The  effects  of  y^* 
the  flow  rate  fraction  of  fresh  feed  stream  have  been  investigated  for 
several  values  of  y^.     The  stability  aiKl  convergence  of  the  computaticai 
became  worse  as  the  value  of  y    became  greater.  Ihis  is  a  result  that  one 
should  expect  fi^am  observation  of  the  boundary  conditions  given  in  equations 
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(87)  through  (90). 

It  should  be  noted  that  in  this  study,  the  values  of  total  holding 
tine,  T]^  and  To  vrere  fixed  cvsn  though  the  fractional  amount  of  fresli  feed, 
"Y^i  ^««s  varying*  To  make  T,  =  8  min,  and  T.  =  10  min.,  either  (aie  can 
specify  the  volume  of  the  second  reactor,  fix  the  flow  rate,  v^-",  and 
adjust  the  size  of  the  first  reactor  and  the  flow  rate,  v^  \   or  one  can 
adjust  the  sizes  of  both  reaotors  and  choose  the  flow  rates  v^  '  and  v^  ' 
according  to  various  cliosen  values  of  y^»    The  results  of  the  computation 
are  i^oun  in  Figures  7  and  8.  The  optimal  solution  and  the  values  of  the 
adjoint  variable  function  for  the  case  y    =  0,5  are  also  shown  in  Table  2 
and  Figure  9  for  the  purpose  of  illustratioi. 
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Table  2. 

Computational  results  for  the  reaction  A-*B-*C  taking  place  in 
two  tubular  reactors  connected  in  series  with  intermediate  feed 

for  the  case  ' 

/2  -  0.5 

Tl 

J3) 

xp) 

43) 

^2 

e(3) 

1.0 

0.15552 

0.67714 

0.00000 

-1.00000 

334,24 

0*9 

0.173^^ 

0.66939 

-0.10252 

-0,93501 

334.^^ 

0.8 

0.19378 

0.65937 

-0.19437 

-0.96982 

334.67 

0.7 

0.21700 

0.64665 

-0.27648 

-0.95436 

334.93 

0.6 

0. 2^^364 

0.63068 

-0.3^^972 

-0.93856 

335.25 

0.5 

0,27^^2 

0.61074 

-0.41484 

-0.92231 

335.62 

0.^ 

0.31033 

0.58588 

-0.47256 

-0.90544 

336.07 

0.3 

0.35270 

0.55477 

-0.523'^ 

-0.38776 

336.63 

0.2 

0.40351 

0.51545 

-0.56814 

-0.86891 

337.35 

0.1 

0.46590 

0.46'+82 

-0.60701 

-0,8^1^34 

338.32 

0.0 

0.54543 

0.39735 

-0.6/+0^;6 

-0,82'^f99 

340.00 

T) 

4" 

4' 

4" 

4" 

e^^) 

1.0 

0.39084 

0.5^*469 

-0.32023 

-0.41249 

327.70 

0.9 

0.41136 

0.52684 

-0.32478 

-0.41045 

327.97 

0.8 

0.43349 

0.50738 

-0.32910 

-0.40834 

328,28 

0.7 

0.45747 

0.48609 

-0.33319 

-0.40613 

328.63 

0,6 

0.48360 

0.46266 

-0.33707 

-0.40383 

329.02 

0.5 

0.51225 

0.43673 

-0.34074 

-0.40140 

329.48 

0.4 

0.54392 

0.40780 

-0.34419 

-0.39880 

330.03 

0.3 

0.57930 

0.37517 

-0.34744 

-0.39601 

330.69 

0.2 

0.61942 

0.33732 

-0.35049 

-0.39294 

331.51 

0.1 

0.66589 

0.29415 

-0.35333 

-0.33943 

332.59 

0.0 

0.72160 

0.24126 

-0.35600 

-0.33540 

334,00 
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k,       COMBINED  USE  OF  TliS  MAXIMUM  P?iINCIPLS  AND  DBIA>3C  PROGRAI-IMING 

The  best  design  job  requires  the  "optttaal"  use  of  every  tool  available* 
kofxng  many  of  the  recently  developed  dynamic  optiiaiaation  techniques,  dynamic 
programming  and  the  raaxiraura  principle  appear  to  be  two  of  the  most  powerful 
for  optimizing  complex  processes  consisting  of  more  than  one  process  unit. 
Each  of  these  methods,  however,  has  comparative  advantages  cmr  the  other, 
depending  upon  the  problem  to  be  solved.  By  oorabining  these  two  techniques, 
it  is  often  possible  to  handle  optimization  of  a  ocraplicated  process  with  less 
effort. 

t»   A  Comparison  of  the  Maxinun  Principle  and  I^mamic  Programming 

The  dynamic  programming  method  reduces  a  eonplicated  problem  into  a 
simple  algorithm  ideally  suited  for  modem  computers.  It  is  based  upon  the 
principle  of  optimality  and  employs  the  technique  of  invajriant  imbedding  [12  j. 
Dynamic  programming,  however,  has  its  limitations.  Since  one  is  in4>edding  a 
i^ole  family  of  processes  in  order  to  obtain  the  solution  of  one  particular 
processf  dynamic  prograraning  is  often  an  expensive  proposition  both  in  time 
and  in  space,  and  since  both  the  computer  time  and  the  required  memory  increase 
rapidly  with  an  increase  in  the  number  of  state  variables,  this  technique  in 
general  can  be  used  successfully  only  for  problems  with  smaJl  nun^ers  of  state 
variables.  !5ven  if  one  has  a  coisiderable  amount  of  computer  time,  the  avail- 
able ra«nory  in  most  modem  computers  still  limits  the  use  of  the  dynamic 
prograzaraing  technique  to  problems  of  not  more  than  three  or  four  state  vari- 
ables. It  is  not  suited  for  solving  continuous  problems  in  i^ich  the  appli- 
cation of  dynamic  programming  leads  to  a  set  of  partial  differential  equations. 
It  is  possible,  of  course,  to  convert  this  set  of  partial  differential  equations 
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Into  a  systen  of  finite  differeno©  equations  arei  apply  the  method  of  dynaiaic 
prograimiiing  to  solve  continuous  problffins.  The  computational  difficulty, 
however,  will  be  evai  more  pronounced  in  this  case  because  a  very  asaall  in- 
creiaent  of  finite  difference  equations  must  be  taka\  in  order  to  obtain  a 
solution  vith  reasonable  accuracy. 

Another  disadvantage  of  dynamic  prc^raming  is  that  it  cannot  be  con- 
v©nl<«jtly  applied  to  processes  in  i^lc^i  the  optimum  caiditicms  at  axiy  stage 
can  be  disturbed  by  conditions  at  the  following  stages.  This  is  a  fairly 
serious  limitation  if  <»ie  considers  the  various  countercurrwit  operations 
employed  in  the  cheiaical  and  petroleum  industries. 

However,  the  definite  and  significant  advantage  of  c^maraio  pr<^rarEning 
in  the  optimlBation  of  processes  tdth  constraints  on  state  variables  cannot 
be  neglected.  Processes  with  bounded  state  variables  give  no  touble  with 
this  method,  because  the  optiiaal  decisicais  deterndrwd  for  the  entire  allowable 
policy  domain  a-atcaaaticaliy  satisiy  the  constraints. 

On  the  other  ha»i,  the  rnaxinwra  principle  can  be  used  to  handl-e  most 
constraints  of  ordinary  optinum  design  jaroblems  with  the  exception  of  those 
in  the  state  variables  which  cannot  be  handled  rout^jwly.  Bf  the  use  of  the 
aaximum  principle  the  optimization  probl<»i  is  essentially  reduced  to  a  maxlHa»- 
or  miniimm-  seddlng  problem  8\ibject  to  a  set  of  ordinary  differential  equations 
and  certain  inequality  constraints.  IMs  set  of  ordinary  differential  equa- 
tions usually  forrss  a  two-point  boundary  value  problan.  Thus  ono   is  still 
limited  by  the  difficijlties  enco^intered  in  solving  a  set  of  simultaneous  non- 
linear differential  equatl<ms  with  two-point  boundary  conditions*  Nevertheless, 
the  nufflber  of  state  variables  that  can  be  treated  b^'-  usually  available  cow- 
poting  equipraent  >^en  the  maxljmim  principle  is  used  may  be  larger  than  that 
vhmx  the  ^mamio  progr^uaning  technique  is  used.  Furthermore,  the  maxiraam 
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principle  is  probably  more  suited  for  solving  c<»itinuous  probleras  than  dynaxaic 
prograj^aair^  is. 

Another  salient  feature  of  dynaraio  prc^rasaning  is  that  it  can  be  used  to 
handle  larocesses  for  which  the  transformation  at  each  stage  is  expressed  in  a 
form  other  than  that  of  a  raatharaatical  equation  (for  example,  a  table  or  a 
graph).  The  raaxiiaum  principle,  on  the  other  hand,  is  applicable  only  to 
processes  with  well-defined  performance  aquations,  and  the  transformation 
functi<ai  must  be  continuously  differentiable  with  respect  to  the  state  vari- 
ables* 

The  COTibined  use  of  the  naodLmum  principle  and  dynamic  prograrandng  can  be 
aniployed  for  the  optimization  of  sosie  complex  processes.  A  process  can  be 
separated  into  different  units,  and  each  unit  can  be  optimized  by  either  of 
the  two  laethods  (or  by  any  other  ajgxropriate  optiEdaatitm  method),  the  over- 
all process  is  then  optimized  by  connecting  the  different  units  by  use  of  the 
principle  of  optlmality.  This  ccMobination  can  be  referred  to  as  the  sequaitial 
tmlom  of  the  maxiauH  principle  and  djmawio  prograoEsir^. 

Ctt  the  other  hand,  since  the  maxiaaua  principle  only  gives  the  necessaiy 
conditions  for  local  raaxiaum  (inininaim),  there  is  no  assurance  that  the  optiiaal 
I^lioy  found  by  the  raaxiHun  principle  is  the  global  optimal  policy.  Sine* 
the  method  of  dynamic  Fnrograraming  emplqjrs  the  so-called  imbedding  technique 
*diieh  is,  in  its  spirit,  similar  to  the  aadiaustive  search,  the  optimal  polity 
obtained  by  dynamic  progranaidng  will  always  be  the  global  optimal  policy 
provided  that  the  interpolati<»i  error  inherent  to  the  method  is  small.  It 
is  thus  proposed  that,  \ih&n  an  accurate  solution  is  desired  or  required,  the 
method  of  dynamic  programming  be  employed  first  to  locate  the  approKimate 
position  of  the  global  raaxiraum  and  that  the  maximum  principle  then  be  applied 
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to  pinpoint  the  maxiMum  point.  This  combination  nay  be  caUed  the  parallel 
union  of  the  maxinma  principle  and  dynamic  programming. 

In  the  ^^ollotfing  section  the  method  of  sequential  union  will  be  presented 
In  detail. 

b.   Sequential  Union  of  the  llaxiinum  Principle  and  Etmasdc  Programming 

•Hie  essential  steps  5ji  the  application  of  the  sequential  union  of  th« 
oaxiBiun  principle  and  dynamic  -rt^rararajjig  are  summarized  as  follows  l21Jj 

(1)  An  entire  process  is  first  divided  into  different  units  (subprocessos). 

(2)  An  appropriate  optirdzation  method  is  applied  to  each  unit. 

(3)  Th«i,  the  optiiffitm  of  the  entire  process  is  obtained  by  the  principle 
of  optimality. 

l^e  separation  of  a  process  into  different  units  is  arbitrary,  and  is 
often  r-ictated  by  a  knowledge  of  the  process  and  its  requirements.  For  example, 
a  chemical  manufacturing  process  idiioh  requires  a  reactor,  an  extractor  and  a 
heat  exchanger  may  be  divided  into  three  or  more  units.  In  the  purification 
of  a  chemical  compound  from  a  mixture  by  several  distillatic«  columns,  one  can 
consider  each  distillation  column  as  a  iinit,  or  one  can  consider  all  the  dis- 
tillation coliimns  together  as  a  single  unit,  depending  on  the  optimization 
requirements  and  on  the  ccraplexity  and  physical  characteristics  of  the  process. 

Each  unit  may  be  ccanposed  of  a  number  of  stages,  or  it  may  Just  represent 
a  continuous  process.  In  a  series  of  stirred  tank  reactors,  each  stirred 
tank  may  be  considered  as  a  stage.  A  tubular  flow  reactor  may  be  considered 
as  a  continuous  process  consisting  of  an  infinite  number  of  stirred  tanks. 

A  two-unit  process  will  illustrate  the  method  of  sequential  union.  The 
arrangements  of  the  units  and  the  number  of  stages  in  each  unit  are  shown  in 
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Figure  10  in  which  Unit  1  is  a  ccaitinuous  process  and  Unit  2  is  a  discrete 
process.  It  is  assumed  that  the  maxicEiia  principle  is  to  be  applied  to  Unit 
1  and  dynamic  prograrming  to  Unit  2,  Ihe  over  all  optiiaua  is  then  obtained 
by  applying  the  principle  of  opbiriality. 

For  the  continuous  simple  process  of  Unit  1,  one  knovs   from  Section  2a 
that  the  perfomance  equations  are  generally  given  by 
dx. 

or  in  vector  form 

^»f(xje)  (91) 

The  aim  is  usually  to  optimize  a  linear  coribinaticai  of  the  final  values  of 
the  state  variables  of  this  unit,  that  isi  to  maximize  (or  minimize)  the 
quantity 

S=  S  c.x.(T)  (92) 

i=a 

with  the  initial  condition  x^Iq)  =  a  given.  To  do  this,  one  generally 

introduces  an  s-dinwisional  adjoint  vector  z(t)  and  a  Plamiltonian  function 

R  which  satisfy  the  following  relations 

s 
H(z.x,t)  =  r  z.S:Ax;Q)  (93) 

i?xl 

and  the  boundary  conditions 

x(tQ)  =  o-  ,   zj^(T)  =  -  c^L  ,    i  =  1 s         (95) 

15ie  optimal  solution  is  obtained  oily  if  the  function  11  is  minimized 
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(or  maxinized)  at  every  monent  t  =  tg  to  t  =  T,  subject  to  aquations  (91) » 
(9^^),  and  (95). 

For  the  discrete  process  of  Unit  2,  the  transforniation  of  the  process 
stream  at  the  nth  stage  mai'  be  described  by  a  set  of  difference  equations  j  7] 

X  "  ^  (x^  •  ^2     »•••#  Xg  t  O-j^t   G2»«»»»  Qy)  t  i  =  1»  2,«»«,  g 
or  in  vector  f orra 

x"  -  ^   (x"-^;  e")  (96) 

^ich  is  the  performance  equation  of  this  unit.  Suppose  that  the  return 
function  (Interval  profit  functiai)  from  stage  n  of  Unit  2  is  h(x:'^  j  6^) 
and  one  wishes  to  raaxinize  the  total  profit  represented  by 

h(x°i  e^)  +  h(x^  e^)  +  ...  +  h(x^-^  e"^)  +  ...  +  h(x^'"^i  q^)     (97) 

It  should  be  noted  that  the  value  of  this  maximun  depaids  only  upon  the 
initial  state  of  tlie  state  variable  and  the  nmdbor  of  stages*    Ihus,  let 

g,,(x°)  =  Max    jh(x°j  0^)+   ...  +  h(x^'"^j  e^^)  I  (98)** 

Using  the  principle  of  optiraality  |  I2J  and  equation  (96),  one  obtains  a 
recurrence  relation  as 

gjj(x^)  »  Max  |h(x°;  e^)  +  g^^^  [f  (x°{  9^)Jj  (99) 

for  N  a  2,  3*«*«*  N«  For  N  =  1,  one  has 

♦  Ihe  superscript,  n,  indicates  the  stage  nuriber. 

♦*  Subscript  !I  indicates  the  total  nuniber  of  stages  but  not  the  Nth 

stage. 
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g3^(x^)  =  Kax  h(x°j  G^)  (100) 

0* 

Thus  the  optimization  problem  of  Unit  2  can  be  solved  recixrsively  by  using 
equations  (99)  and  (100). 

In  order  to  obtain  the  over-all  optiimim  for  the  process,  the  optinua 
of  Unit  2  naisFt  be  achieved  first  by  the  method  of  dynardc  prograwaing  1 12  U 
Since  Unit  1  is  to  be  optimized  by  the  naximuin  principle,  special  attention 
must  be  given  to  the  objective  function  for  this  unit.  Because  caie  has 
alreacfy  obtained  the  optimum  of  Unit  2  by  dynacdc  pxrc^ranraing,  equation  (92) 
can  be  modified  to  include  the  effects  of  Unit  2.  It  should  be  noted  that 
the  optinum  retxim  for  Unit  2  is  in  goieral  given  in  tabular  form.  For 
every  possible  conbination  of  the  initial  condition  x°,  the  value  of  the 
optimum  return,  g^^  [^see  equation  (99) j,  is  listed.  In  order  to  use  the 
function  g.j(x  )  in  establishing  the  return  function  for  Unit  1,  it  must 
first  be  expressed  in  analytical  form.  This  can  be  done  sometimes  by  using 
appropriate  correlation  methods.  For  example,  a  polynomial  approximaticai 
appears  to  be  promising  in  these  correlaticais.  In  a  one-state  variable 
process  the  gjj(x°)  may  be  represented  in  an  analytical  form. 

Vhen  the  function  gfj(x°)  is  expressed  analytically,  the  objective 
function  for  the  entire  process  can  be  obtained.  If  one  lets  f  h(xj  8)dt 
be  the  return  (profit  or  cost)  from  Unit  1,  the  objective  fimction  for  the 
overfall  sj'-stcan  is 

S  =  gjj  j_x(T)J  +  r  h(x!e)dt  (101) 

xtere  x(T)  -  x  ,  that  is,  x(T)  represaits  the  final  value  of  the  state  vector 
of  Unit  1  and  is  equal  to  the  initial  value  of  the  state  vector  of  Unit  2, 


Equation  (101),  v^iich  corresponds  to  the  recurrence  relation,  equation  (99) » 
is  established  by  the  principle  of  optimality.  '£he  method  outlined  by 
equations  (91)  through  (95)  t  niay  now  be  modified  to  take  into  account  the 
objective  function,  equation  (101).  To  accomplish  this,  a  new  state  variable 
is  introduced  as 

T 

or 

Vl^^^"  ^L^'^^U+I  h(xje)dt  (103) 

and 

The  problem  is  now  transformed  into  that  of  maximizing  the  final  value  of 

X  !•  The  Ilamiltonian  function  H  in  eqxiation  (93)  is  modified  as 

s+1 
H(x,2,e)  =  E  2ifi(x;0)  (105) 

i»l 

and  the  new  adjoint  variable  z  ,  is  introduced  as 

%i=.^=0  (106) 

The  boundary  ctmditions  for  the  adjoint  variables,  equation  (105),  should 
be  changed  to 

Z|^(T)  =  0  ,        i  =  1,  2,,,,,  s 


(107) 


23,3_(T)  =  -  1 


Ftoti  equations  (106)  and  (10?),  one  obtains  8-^.1  (^)  *•  -  If  and  hence  one 
has 
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s  ^1 

n(x,a,e)  =  E  z^f^ixie)  -  -gfii  (108) 

All  other  relatlwships  of  equaticaig  (91)  and  (9^+)  ranain  the  sane.  This 
optimization  problem  can  now  be  solved  by  using  methods  introduced  in  Section 
2  for  the  case  of  a  simple  process,  i.e.,  one  involving  no  separating  and 
combining  points.  Cttce  the  optimization  calculaticm  for  Unit  1  is  completed, 
the  final  values  x(T)  are  obtained.  These  values  provide  the  optimal  initial 
values  for  Unit  2,  The  optimal  policy  for  Unit  2  can  be  recovered  by  using 
equatiws  (99)  and  (100). 

So  far  a  system  vriiose  first  unit  is  optimized  by  employing  the  maxiniura 
rrinciple  has  been  dealt  with,  '^en  the  first  unit  in  the  sequence  is  not  to 
be  optimized  by  the  maximum  princii>le,  the  same  procedure  is  still  applicable. 
However,  one  must,  in  general,  calculate  a  grid  of  values  associated  ^d.th  the 
unit  optimized  by  the  maximum  principle,  since  in  this  case  the  initial 
conditions  of  the  units  using  the  maximum  principle  are  not  given.  Generally 
there  are  no  givwi  initial  conditiais  for  any  unit  except  the  first  one. 
''Sien  the  grid  of  value  is  obtained,  the  optimal  policy  of  the  units  preceding 
this  unit  can  be  obtained  by  applying  the  principle  of  optimality.  Note  that 
if  the  unit  which  is  to  be  optimized  by  using  the  maximum  principle  is  a  one- 
state  variable  process  and  the  functional  dependence  of  the  optimal  value 
vith  respect  to  the  feed  conriitiOTi  to  this  unit  can  be  6btain<xi  analytically, 
the  calculation  of  a  n^id  of  values  is  not  necessary. 

In  ^at  follows,  t^'W  examples  will  be  presented  to  illustrate  the  use 
of  the  method  of  sequential  unicm.  The  examples  are  chosen  so  that  the  over- 
all optimum  can  be  obtained  from  the  analytical  solution  without  any  calcula- 
tion of  a  grid  of  values. 
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«•   Tab'jlar  Flow  Reactor  Followed  by  a  Multi-ertage  Cross-currwit  Sxtraction 
Process 

In  this  section  a  two-unit  process  i^ich  consists  of  a  tubular  reactor 
of  length  IL  (continuous  process)  followed  hy  a  two-stage  cross-current 
sxtraotor  (discrete  process),  as  shown  in  Figure  11  la  considered. 

It  is  assumed  that  an  autocatalytlc  reaction  of  the  type 

B   R.^R   R 

takes  place  in  the  tubular  reactor.  TUxe  product  stream  is  then  fed  to  the 
liquid-liquid  extractor  for  purification.  It  is  assumed  also  that  reactant 
B  la  injdseible  in  the  solvent  and  only  the  product  R  is  extracted*  It  is 
desired  to  maximize  the  over-all  net  profit  by  finding  the  optimal  t«aperature 
profile  along  the  txibular  reactor  and  by  allocating  the  optimal  amount  of 
mailing  irater  at  eac^  of  the  tvro  extractors. 

Unit  2  will  be  optimized  fiirst  hy  using  dynamic  programming  and  th«» 
Unit  1  by  the  maximum  prlncipCLe.  The  over-all  optimiara  is  then  obtained  by 
use  of  the  principle  of  optimality.  Die  rate  of  the  autocatalytie  reaction 
can  be  writtwi,  with  respect  to  production  of  reactant  R,  as 

where  G  is  the  molar  concentration  of  componwt  ij  b  and  r  are  the  order 
of  reaction  with  respect  to  components  B  and  R  respectively  in  the  forward 
reaction  and  s  is  the  order  of  i*eaction  of  the  reverse  reaction.  ?>rom  the 
stoichiometry  one  has  at  evexy  point  alor^  the  reactor 

where  C  .  is  the  initial  molar  ocxicentration  of  reactant  i.  "nierefore, 
equation  (109)  can  be  written  as 
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Rftarranglng  the  above  equation,  cne  obt&ins 

Bte  stead^'-state  differential  laaterial  balance  for  reactant  R  in  the 
differwitial  section  di*  of  reactor  yields 

Rate  of  input  of  coraponfflit  R  at  the  inlet  of 

{  i 

^dif fer«itial  section  de  of  the  reactor  line 
Rate  of  output  of  componcmt  R  at  the  outlet  of 

Mifferantial  section  dl  of  the  reactor  line   ^ 
Rate  of  production  of  coojponent  R  in  the 

Mifferential  section  di  of  tee  reactor  Uxor 

m   0 


i*e* 


W-K  -^^i  ^H^}-" 


_^  -  ^  W 

vhere  v  is  the  volumetric  flot?  rate  in  the  reactor  and  A  is  the  cross- 

sectiwial  area  of  the  reactor. 

Substituting  equation  (HO)  into  this  equation  yields  «*..,« 

^  -  7  L^i^^  •  y  §)'  CCo)^^'.  k2§)"  io,fj  cm) 

^i 

In  terms  of  fractional  molar  concentration,  Xj^  =  ^»  <»e  can  rewrite  equation 
(lU)  as 


5S 


^ = A  L'',  (1  -  .^)\x,nc,)^^  -  ^,(^f(c,f-'-j     (U2) 

liiere  subscript  1  refers  to  reactant  R  (ccmiponent  !)•  If  one  lets 
Cq  =  Cg  +  Cfj  =  1,  equation  (112)  can  be  replaced  by 

^  =  A  Lki  (1  -  ^)\x^f  -  k2fei)=J  (113) 

Ihis  is  the  parformance  equation  for  this  unit. 

One  derives  the  transformation  equaticai  of  Unit  2  by  considering  it  as 
an  n-stage  process  (see  Figure  12),  The  material  balance  i-dth  respect  to 
product  R  around  the  nth  stage  gives  I  22J 

wV  =  q(xj-^  -  Xj)  (114) 

\iiBre  V  ~  wash  water  ratot  Ib/hr 

y  =  coi^x>sition  of  product  R  at  Hie  outlet  of  the  extractor,  dimensionless 
q  s  vp  ,  mass  flow  rate,  Ih/hr 


or 


p  =  density  of  mixture  (assmned  to  be  constant),  — » 


n  n   n-1   n 
u  y  =  :i^   -  x£ 


(115) 


where 


u  »  W  /q  (116) 

Let  Q„  =    2    ypy'^  =  q(x5  -  s^)  be  the  total  aiaount  of  R  extracted  and 

N      ^^ 

Gj.  »  ^  -^  be  the  total  amount  of  washing  water  supplied.  Ihen  the  net 

profit  function  for  Unit  2  nay  be  i^rritten  as 

P„  =  %  -  ^Gj,  (117) 

'twere  X  is  the  relative  cost  of  washing  water  and  it  can  be  regarded  as  a 
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Lagrange  multiplier.  Squation  (117)  is  equivalent  to 

R,  -  q  I     (x?-^  -  x?)  -  X  I    ^  (118) 

^^    n=l  n=l 


(119) 


Substitution  of  equation  (116)  Into  equation  (118)  gives 

tfcich»  together  vith  equation  (115)  are  the  performance  equations  of  Unit  2, 

Now  consider  the  return  function  of  Unit  2.  Frm  equation  (119) »  the 
return  function  of  Unit  2  is 

f„(x?)=  Max  I    u"(y"-X)  .    4  =  *  ^^°^ 

^      |u"|  "=^ 

If  a  linear  equilibrium  relationship  is  assumed  in  each  stage,  the  relation 
between  the  ccmcentrations  x  and  jr  is 

X  =  oy  (121) 

Substituting  equation  (121)  into  equaticai  (120)  yields 


fjj(a)  =  Max  E  u"  (^  - 


^^nj.n=l 


\)  (122) 


According  to  equation  (99),  one  has  the  following  recursive  relation  as 

fy(a)  =  I'lax  {v?-  if'\)+  ^.iCxJ))  (123) 


ul 


with  :^  =  a  and  fQ(a)  =  0, 

For  a  tvio-stage  extractic«i  process  (K  =  2),  cor^ining  equations  (115) 
and  (121),  gives  for  the  second  stage 

:C  =  xT  -  u  V  =  07^ 


5? 

or 

U+Of 

(12i^) 

Bquation  (123)  f or  N  =  2  is 

1  ^ 
fgCa)  =  Max  \n^  i-f  - 

X)  +  fj,  (xj)j 

(125) 

Tdiere 

2 

f,(4)  =  ^^ax  tu^  (?  - 

u 

X)  ^  fo(xJ)) 

(126) 

Substituting  oquaticai  (124)  Into 

equation  (126) 

and  noting  that  £q(,^)   =  0, 

one  obtains 

f,  (J;)  =  Max  iv?   (-S- 

u 

-uj 

(127) 

Partially  difforaitiating  the  quantity  in  the  bracket  of  equation  (12?) 

2 
with  respect  to  u  and  setting  it  equal  to  zero 

"7 
to  solve  for  u  , 

one  obtains 

1      [4°' 

u  =J  -p-* 

(128)» 

Substitutinc  equation  (128)  into  equation  (124), 

one  obtains 

'^^//^ 

(129) 

Substituting  equation  (128)  into 

equation  (12?) 

fax 

gives 

(130) 

fl(4)  =^  -zj  c^X 

*     Ihe  upper  bar  indicates  the  optteal  solution. 
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provided  that  u^  >  X«,  This  is  the  optimal  interval  return  function  for 
the  second  stage  of  the  extracti<»i  process  with  respect  to  the  outlet 
conditicai  of  the  first  stage.  Similarly,  for  the  first  sta^e,  applying 
equation  (125).  one  has 

fgCa)  =  Max  \^   (9^  -  3^)  +  \^^)\  ^^1) 

Substitutinc  equation  (130)  into  equation  (131)  yields 

f2(a)  =  Vox  \^   (~  -  X)  -J-  3^  -  2/°^^  ^ "'^l  ^^^^ 

Corabini):^  equations  (115)  and  (121),  for  n  =  1,  one  has  for  the  first  stage 

x£  =  x^  -  uV  =  oy 
or 

u  +  0-   u  +  a 

Substituting  equation  (133)  into  equati<ai  (132)  and  si^-nplifyinc  the  resulting 
equation,  one  obtains 

f2(a)  =  Max  {a  -  tt\  -  2   /i|i-^  +  a\j  (13^) 

^1  W    U  Q' 

Partially  difforentiatir^  the  quantity  in  the  bracket  of  equation  (134) 
with  respect  to  u  and  setting  it  equal  to  zero  to  solve  for  u  ,  one  obtains 

1 

2  5 

?-(-^)    -(^  (135) 

Substituting  equation  (135)  into  equation  (133)  i/delds 


a 


Substituting  equations  (135)  and  (1}6)  into  aiuation  (13^),  one  has 

1 

5 

fgCa)  =  a  -  3  l^^fi^f&j       2ctX  (137) 

Substituting  equation  (136)  into  equation  (128)  and  noting  eqtiation  (135) i 
one  obtains  the  optical  policy  for  the  two-stage  extraction  unit  as 

1 

—     2  5 
?  =  vi2  =  (S2IA)     -  a  (138) 

TJiis  neans  that  the  wash  xrater  is  divided  equally  between  the  two  stages* 

To  obtain  the  optiroun  for  the  over-all  process  or  the  niaxiinuia  coribined 
profit  fron  the  reactor  and  purification  sections  t  the  objective  function 
must  be  established  for  the  entire  process.  Since  no  product  R  can  be  sold 
until  it  is  porifiod,  there  is  no  direct  return  available  frora  Unit  1.  The 
operating  cost  for  Unit  1  is  assumed  to  be  poroporticnal  to  the  amount  of 
reactant  B  converted  into  the  product  R  as 

P.(xj^(L)  -  Xj^(O))  (139) 

^ere  3  is  a  proportionality  constant  and  x,  (L)  and  x,  (0)  are  the  fracticaial 
mass  concentrations  of  R  at  the  outlet  and  inlet  of  the  reactor  respectively. 
Since  the  raaxiimun  profit  frcsn  the  purification  section  has  already  been 
established,  the  over-all  net  profit  function  can  be  written  frot;  equations 
(137)  and  (139)  accordinc  to  the  principle  of  optiioality  as  I  see  equation  (lOl)J 


S  -  4  -  3[(<v)^(X)2xJj  "+  2C.X  -  0J_x^(L)  -  x^(0)J         (140) 
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Remerabering  that  Xj^^W   =  a^t  o^®  can  revorite  equation  (l^i-O)  as 

S  =  Xj^a)  -  \^icxf(\)\iL)^  +  2aX  -  -^^Xj^iL)   -  x^iO)^ 

(m) 

To  find  the  docisi<ai  functicai  OU)  >diich  aaxinizes  the  objoctive  function. 

on©  first  introduces  a  second  state  variable  x^  defined  as 

1 

5 

(1^2) 

Difforaitiation  of  ©quaticm  (1^1-2)  with  respect  to  £.  gives 

^  =  (i-P)  ^  -  (o.x)5  >^  "5  *^ 

(l'+3) 

This  equation  and  equation  (113)  fora  the  ccxaplete  set  of  performance  equations 

for  Unit  1,  The  initial  conditions  are 

x^CO)  =  Y 

a^) 

1 
XjCO)  =  Y  -  3  L(a)^(X)\j  +  2a\ 

(1^5) 

Substituting  equation  (113)  into  equation  (143),  one  obtains 

2    2 
^  =  [l^(c.X)^  T^  ^Jk^a^Xjf{x^f^k^(x^f^ 

(1^) 

According  to  equations  (7)  and  (8),  we  can  -write  the  ilaniltonian  H 

and  the  adjoint  variables  for  this  imit  as 

H  «  -^  j_kl(l-.XL)^(x^)^-k2(x^)Sj  -^  ^  L^(l.X;^)^(x^)^-k2(x^)Sj 

2^2 

(1^7) 
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5ji  =  -  -^  LblCj^(l-X;^)^^(x;^)''+ki(l-X3^)^rCx^)^^-k2s(x;L)^\ 

.  ^2  j_bk3_(l.Xi)^-^(xi)'>k3_(l-Xi)^(xi)^^-k2s(x3_)«-^J 

2    ^2 

-  ^      - 

im) 

=  0 

iV^9) 

Sine©  the  objective  ftmcticm  is  XgCL),  xw  have  On  =  0  and  Cg  =  1.  According 

to  equation  (95) t  one  has 

z^iL)  =  0 

(150) 

z^(L)   =  -  1 

(151) 

Corablnine  eqtiation  (151)  viith  equation  (1^9)  yields 

ZgU)  =  -  1 

(152) 

Substituting  equation  (152)  into  equation  (14?) i  one  has 

2    _  2 
H  =  1  Lk^(l-X3^)^(x^)''-k2(x3_)^J  {z^^  ll^^iaXp   (x^)  ^_[| 

(153) 

Vhen  the  Blnlimini  of  H  occurs  at  the  stationary  point,  one  can  obtain  the 

optimal  decision  functicm  9 (A)  by  partially  differentiating  equation  (153) 

vith  respect  to  e  and  setting  it  eqiial  to  aero,  that  is 

,^^.i^.l.^.J^.,ih^,,,^,.,^,r^.,^,s^^ 

=  0 

dkn    dk^ 
T,iiere  ^«r=  and  -w^  can  be  obtained  by  differontiatinc  the  equations 

€k 

-     2 

with  respect  to  9.     One  thus  finally  obtnins  the  optimal  ten^serature 

.  o(i). 

«8 

Tf,\         f"^        ~'?s                         }■                       , 

(15^) 

^"'-*      K         ta  LE(Kir^(l-Xi)*J 

\here 

_      ^0  ^ 

^"^0^ 

Substituting  equatiwi  (15^)  into  equation  (113)  yields 

ar "  ^  i%o(i-»i>^-i'""L5(==i>'"'(i-^i)''j  ^ 

-  ^.'So  L'^(=i)'"°G-==i)''/^) 

(155) 

lAiere 

1      Eg-a^                           2      ^-% 

In  integrating  equation  (155)  the  activatioii  ^oi^ies  arc  ahosen 

in  order 

that 

■^  =  2          or          ^1  =  ^          ^^          ^2  '^  ^ 

(156) 

It  should  be  noted  that  oniiation  (155)  is  valid  for  any  value  of 

1 

and  Eg. 

Equation  (156)  is  choson  for  t}io  stniplificaticai  of  the  probleia. 

Substituting  aquation  (I56)  into  equation  (155)  and  simplifying  the 
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resulting  equation,  one  obtains 

^  =■  I  L"-=10'  -  '^o(?)'j  Ltt-=^)*(--=i)'""=j  (157) 

If  one  assumes  that  the  reaction  is  first  order  with  respect  to  each 
con^nent  for  the  f onrard  reaction  and  second  onler  with  respect  to  the 
reversible  reaction,  that  is,  b  =  r  =  1  and  s  =  2,  equation  (157)  can  be 
Integrated  as 

The  integration  constant  C  is  evaluated  by  using  equation  (1^)  as 

«i-  =  C  (159) 

1-Y 

Substituting  equation  (159)  into  eqioation  (158)  yields 

I  =      .,     '^    .       ■   (1     -  -L)   •.  (160) 

Con&ining  equations  (15^+)  and  (160)  one  obtains  the  optl'ial  mlicy  for 
Unit  1. 

Solving  equation  (l6o)  for  x,(L),  one  has 

or 

Iy+  i^dc^  -k  f)(i-y) 
X.  (L)  =  ^ -^ — 22 (161) 

Substituting  equation  (I6I)  into  equations  (139)  and  (137)  gives  the  optlioal 
profit  (or  cost)  fear  the  reactor  and  the  purification  section  respectively, 
Ihe  over-all  optiinal  net  profit  is  obtained  by  substituting  equation  (I6I) 
into  equation  (141). 
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d.   Continuous  Stirred  Tank  -leactor  Follc.^  by  a  Tubular  Flow  Reactor 
with  Product  Recyclo  to  an  liitemal  Point 

In  this  sGOtion  a  chonical  procecn  consisting  of  a  continuous  stirred 
tank  reactor  folloi-red  by  a  t'lbt-xlar  roaotor  :olth  product  r'^cycle  to  an  inter- 
nal point,  as  shown  in  Figure  13 »  \A2X  be  considered.  Unit  2  (the  tubular 
reactor)  consists  of  four  branches,  tiio  food  line  of  lenjtJi  L.  (brandi  1), 
tho  ror.otor  of  leiicth  L^  Coranc]-!  2),  the  diGc!iar:;e  line  of  longtii  L,  (branch 
3)t  and  the  rec^rcle  line  of  length  I^^  (branch  '0  vdth  voluraetric  flow  rates 
of  v^-'-S  v^^',  v'^\  and  v^^^  respectively.  Tn  Unit  1  (CSTR),  the  volunetric 
flow  rate  and  the  volune  of  reactor  are  given  as  v  and  V  rospectivoly. 

A  first-order  reversible  chemical  reacticai 


^3 


is  taking  place  in  this  series  of  reactors.  It  is  desii^  to  inaxirnize  the 
conversion  of  product  B  at  the  end  of  tlie  tubular  reactor  by  the  proper 
nanipulatiOTi  of  the  temperature  in  the  stirred  tank  reactor  and  the  temperature 
profile  along  the  tubular  reactor. 

First,  the  luaximun  principle  is  applied  to  optimize  the  tubular  reactor 
unit  and  the  over-all  optiriun  is  then  obtained  by  the  principle  of  optiriality. 

Let  us  establish  the  perfomance  equations  for  both  units  as  follows; 

A  steady-state  rAterial  balance  for  the  ith  ccnponent  over  a  Z"TX  is 
given  by 

■irate  of  input  of  the  ith  c<xipQnent|  +  irate  of  production  of 
the  ith  conponent|  =  -jrate  of  output  of  the  ith  ccaiponent|- 
that  is. 
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VCJ+R^V  = 

< 

(162)« 

•iiheve 

C.  =  inlot  Tiase 

c-onoontration  of  the  itla  coinponeat, 

lb 
ft3 

C^  =  outlet  mass 
V    =  total  volun 

!  conceiitration  of  1th  coEpanant      , 
letric  flow  rate                                   , 

lb 

■ft3 

ft3 
hp 

"    "  nans  react! 

X 

.on  rate  of  i)ro<&iction  for  itl> 

conp-taiont 

t 

lb 

ft^-}ir 

V    =  volume  of  reactor                                                   , 

ft3 

For  the  reaction 

^1 

^2 

in  \fhAch  both  the  fonrard  and  bnckward  reaction  are  assumed  to  be  first- 

order,  the  reaction  rate  R^  in  terns  of  the  production  o^ 

coi.ponont 

2  is 

Given  by 

h  =  h'^A ' 

.  kC^ 

(163) 

\iiere  the  specific  reaction  rate  k.  foUoira  Arrhenius  latj 

r  as 

^1  "  ^10  ® 

(164) 

lAiere  the  reaction  tei^peratiu^  9     is  the  decision  variable.     All  other 

*    Superscripts  indicate  the  stace  in  Unit  1. 
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notaticai  is  defined  in  Soction  3a.  Substitution  of  aqr-ation  (I63)  into 
equation  (l62)  for  conponent  3  Gives 

YC^  +  (k^c^  -  k^cj)  V  =  V^ 

Dividing  the  above  equation  by  v(Ct      C~),  one  obtains 


Assumlmc  that  Gt  C_  =  C~-f-Cris  constant,  and  letting  ^  =  t,  the  preceding 


A  '      B        A        B 
equaticMi  becones 


y^+  tik^^  ^  k^)  ^  :^  (165) 


v^ere 


0  4 

4    ^ 

1    i 

■^ '  o'A ' 

^+4 

V°B 

In  aquation  (I65)  :Xri  and  Xg  ropresoiit  tie  fractional  mass  concentration  of 
3  and  A  rospoctivoly  at  -tJ-io  outlet  of  CjTR  and  t  is  tiie  holding  time  in 
the  reactor. 

SincG  ::C  +  x^  =  1,  equation  (I65)  can  be  replaced  by 

T     x^  +  k^  (0^)  t 

4  — ^^ — h r-  (^^) 

■"  1  +  t^k^Co  )+k2(e-')j 

Fx^ation  (166)  is  the  pci-formanco  oq-.'.ation  for  Unit  1. 

A  steady  state  material  balance  Tjdth  respect  to  couiponont  3  over  a 
differential  section  d,?  in  brancii  2  of  Unit  2  leads  to  (see  Section  ^c) 

^^^  V'-^  (2) 

V   +  V 
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vtiOTB  A^^'  =  cross  sectional  area  of  tubular  reactor  (Branch  2) 

v'^'  a  voltuaetric  n.ow  rate  of  feed  lino 

v^^^  =  voluraetric  flow  rate  of  recycle  line 

Since  no  transformation  is  taking  place  and  no  control  action  is 

tmpoiwd  on  branches  1,  3,  and  4,  one  has  the  following  perfoiwance  equations 

d^W   dx,»>   drW 

(168) 

for  the  feed,  discharge  and  recycle  lines,  respectively. 

Bie  initial  state  of  Unit  2  ia»  according  to  toe  above  equatiai 

x^^5(o)«x^(l)(l^)-x^l 

(169) 

Conservati<ai  of  mass  at  the  coaabining  point  requires  the  junotioai  equatiaa 

to  be 

,^'^4'ho)  .  v»^"(ia)  *  v<'*^'*>(i^) 

or 

x{25(0)=<«^&^Il)+tef»>(V 

(170) 

^ere 

^(1)              ^(^) 

"  '  ^(x\^Lk)  '             ^  "  ^(l)^^^> 

The  jimctioi  equatirais  at  the  separating  point  are 

(171) 

By  substituting  equations  (163),  (169).  and  (171)  l«to  equation  (170) 

one  Obtains 

x[^^(o)«c«5;f  ^x^^hi^) 

(172) 
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It  Ifi  desired  to  aaxlmize  the  yield  of  product  B,  that  is,  to  maximize  tii* 
objective  function 

S  -  xp^I^)  -  4 
Bie  dbjeetive  function  in  equation  (173)  can  be  written  as 

\iaer9 

Che  can  ccnsidwr  that  S-  and  Sg  are  the  Interval  porofits  trcm  Unit  1  and 
Unit  2  respectively.  Thus  the  ppdtilesa.  becomes  one  of  maximizing  the  total 
profit  (S,  +  S,)  by  a  proper  choice  of  the  toiqwrature  in  the  CSTR  and  the 
tubular  reactor.  According  to  equation  (99)  •  one  can  write  the  recursive 
relation  as 

S^(7^)  =  Max  LaL(4;9^)  +  g^(4)J  (176) 


}ii9iPB 


g^CxJ).   Hax  LS2(3^ie(2)(^jjj„   j^  L'i^^CV-'ij   (177) 


Unit  2  will  be  optimized  first  by  using  the  maxinura  principle  and  the 
overxall  optimum  can  then  be  obtained  by  the  use  of  equation  (176)«  la 
optimizing  Unit  2,  one  ctmsiders  the  feed  c<mcentratian  to  the  tubular 
reactor  to  be  givant  for  the  time  being  i  that  is»  s^  «  a  =  constant. 


n 


Thus,  naxlraiaing  the  interval  profit  represented  in  equation  (175)  Is 
equivalent  to  laaxlBizlng  x^Hl^)*    ?y  "Uiis  consideration,  the  optiMzation 
of  Unit  2  reduces  to  the  problem  of  a  sic^xLe  recycle  process. 
Accordine  to  equations  (7)  and  (8),  aie  has 

H<«  -  h'^'  .  H<*'  .  0  (178) 


6^^      d^3)      ^^) 

ir""'sr' '"UT 


(aBo) 


dasi^^  »(2)  ^c^ 

V     +v 

Since  the  objective  function  is  x^  0^)»   one  has  c£  '  »  1.  Ihus,  at  the 
final  point » 

4^\l^)  =  .  1  (182) 

(2) 
If  one  assumes  that  the  ndninam  of  H  occurs  at  the  stationary  point  (6^  '(i) 

is  not  constrained),  one  can  detendne  the  optimal  ternperature  profile  by 

taking  the  partial  derivative  of  H  with  respect  to  6  and  setting  it  equal 

to  zero.  Thus,  fran  equatiai  (179),  one  obtains 

^(2)     ^(Z\^)  ak^    dk2    (2)   *1  ,   „   ,,o,, 

;^=- 7177757  L<;;T5T  +  y!T''i  -^J-"  <^3) 


vhere 


a  -  (1») 
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solving  equation  (183)  for  6^^^  yields 

Q^^hl)" ^'%)  C185) 

vhere 

.      ^20^ 


ho  h 


Substitution  of  equation  (185)  into  equaticn  (I67)  yields 
dx{^^      .(2) 

^diere 


(186) 


In  the  integraticai  of  equation  (186),  the  activaticsi  energies  of  the  forward 
and  badcwaird  reactions  have  been  chosen  such  that 

^  «  2 

or 


Substituting  theso  values  of  X^  Into  equation  (186),  and  perforrdng  the 
Integration,  one  obtains 

a)    (^) 


V  tv 


A^^^dcj^r^.k^o?-")'^-^! 


|_--2^  +  ln(l-xpbj  -  A   C      (189) 


Cne  evaluates  the  integration  constant  C  1^  the  use  of  the  initial  caidition 
given  in  equatloi  (172).  Thus, 

Substitution  of  the  above  equaticsi  into  equation  (189)  yields 
_(1),_(^)      .1         (2)  X 

-  3n  [l.(<»^^px[^^  (I^) ) j}  =  je         (191) 

Equation  (191)  indicates  the  optimal  concentration  along  the  tubular  reactor. 
At  the  end  of  the  reactor,  that  is,  at  i  «  l^t   equation  (191)  Involves  oily 
one  unknoxtfi,  x[^^(I*j).  It  Is  seen  that  equation  (191)  cannot  be  solved 
eoqplicitly  for  x[^'(I^)*  Di  other  words,  one  cannot  express  gj^Cn^l)  [^see 
equation  (177) J  explicitly  as  a  function  of  x^Hl^)  or  35J.  Ihe  optiraua 
value  of  x^hl^)  with  regard  to  xj,  the  inlet  condition  of  Unit  2,  la 
implicitly  given  b7  the  following  equation 

C«  {-i*  t  In(l-y) \~^  -  In  Ll-(a2jtPy)j}  -  I2  »  0     0-92) 

^"■y         l«»(aacT+  Py) 

m 

F  (xj.y)    =    0  (193) 

■rfiere  y  represents  the  optinam  value  of  X£  '  (Lj)  and 


n 


HencOf  one  finally  has 

For  the  pptisdssation  of  Unit  2,  3C  \ms  considered  as  a  given  initial  condition, 
the  functional  interdependence  &,  (xJ)  is  ifflfplicitly  given  by  equation  (192) 
or  (193).  Substituting  equation s(17<f)  and  (19^)  into  equation  (I'A).  aw 
obtains 

g2(:^)  »  Max  l4  -  4 -^  y.3tJj=MaxL7-3^  OSS) 

The  quantity  in  the  bracket  is  a  function  of  6  and  x^,  because  y  is  a 
function  of  jdj*  tdiichf  in  turn,  is  a  function  of  9^  as  represented  in  equation 
(166)*  These  functional  relationships  idll  be  denoted  as 

GiT^lQh  =  y  -  xj 

Wben  the  taaxiPtam  value  of  0(xPte^)  occurs  at  the  stationary  point  with 
respect  to  8  ,  the  optical  value  of  9  can  be  obtained  by  setting 

^  «  dy^  .  ^  „  0  (196) 

ae  63^    de 

Cfte  evaluates  ^s*  by 


-^|r-  (197) 

According  to  equation  (197)  t  Qu©  obtains  trom  equation  (192) 


1C 


^"  (l-y)\l.(o«J+Py)/ 

Substituting  the  preceding  two  eqviations  into  oq\iation  (197)  t  <»o  obtains 

-  cr(l.f)2(a5A.3f) 

SL. « 5 ^      g     H (198) 

(bj      yil-(axit3y)j    -  P(l-y)  («>^  +f3y) 

Differaatiating  equation  066)  with  respect  to  G"^  yleMs 


—i  t  f  (t)2  (ko  -4 


dk. 


1    -4 1  f  (t)2  (kg  -4  -  k,  -4) 


"7[ 


Ll+(ki+k2)tJ^ 


Rflcalling  tliat  the  faacticais  k.  (0  )  and  ^i^  )  are  given  in  eqtiation  (iS^f), 
the  differentiation  of  the  pirecedlng  equation  would  becor» 

d^      ^LV      ^'ttk3,(V^)e        ^^^'j 
d?°  I  l+k^(e^)  -^  k2(e^)t  p 

Substituting  equations  (198)  and  (199)  into  equation  (196)  yields 

ae^        R(e-'-)\i+(k3;fk2)t/ 

<y(l-y)^(o4^Py) 

y  Ll-(caJ^-9y)/  -  &(l-y)^(QXi+3y) 
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Since  oaC+^f  4  0  and  1-y  4^  0  for  any  conversion  less  than  100  per  cent, 
the  only  possibility  for  eqxiation  (200)  to  be  equal  to  zero  is 


6^ 


When  the  value  of  teiaperattire  must  be  finite,  the  latter  case  would  give 
the  optimal  temperature  as 

e-*-  a £-^ (201) 

Rln  (^) 

where  t  is  the  holding  time  in  Unit  1  defined  as  s  and  X<t  »  w  ■■A  » 

It  can  be  shown  that  equation  (201)  gives  the  saioe  optljoal  teiaperaturt 
as  one  will  have  for  a  single  CSTR.  Thus,  the  optimal  docisicai  to  be  made 
at  the  CSTR  r«nains  the  sane  or,  in  other  words,  it  is  independent  of  the 
downstream  condition.  This  is  actually  observed  fnaa  equation  (200)* 
Equating  equation  (200)  to  aero  is  equivalent  to  equating  equatiai  (199) 
to  zero,   ^i<A  yields  the  optimal  decisicai  for  a  single  CSTR« 

C*ie  comes  to  the  conclusion  that  the  optimal  policy  for  this  ti«>-unit 
^stem  is  to  choose  the  tengjerature  which  will  raaxirdse  the  outlet  concen- 
tration of  Unit  1,  and  let  this  outle'  concentraticn  be  the  given  inlet 
conditicai  of  Unit  2  and  then  maximize  the  outlet  concentration  of  Unit  2  by 
the  proper  choice  of  tengieratur©  al<»ig  the  tubular  reactor.  These  facts  may 


also  be  verified  by  the  following  reasoning.  H^uation  (198)  indicates  how 

iima] 
.(2) 


x^  affects  the  optimal  final  concoitratioQ  f.    Equation  (189)  indicates  that 


100^  conversion  (x^  '  «  1)  gives  rise  to  a  reactor  with  infinite  length. 


n 


that  is*  the  oawersion  laaat  be  less  than  100^.    Therefore,  the  fc^owing 
inequality  rauat  hold  for  the  numerator  of  equation  (198) 

oa-y)^(axJ-f  Py)  >  0  (202) 

Furthemore,  m  have  y  >  (oaJ+  0y)  and  p  <  1|  hence,  the  following  inequality 
fflist  hold  for  the  daiocdxiator  of  equaticw  (198) 

y  Ll.(caJ+  ey)  J  >  P(l-y)^(aa^+  Py)  (203) 

Since  both  the  numerator  and  dencMdnator  of  equation  (198)  are  positire,  one 
has  the  following  Inequality 

>  0 

33ti  other  words,  y  is  an  increasing  function  with  respect  to  x£.  In  order 
to  obtain  a  hi^er  value  of  y,  one  must  have  a  hi^w  value  of  xj  regardless 
of  the  operating  policy  in  Unit  2«  this  simple  but  surprising  result  appeara 
to  be  a  consequence  of  the  characteristics  of  the  reversible  reaction  A  ^=^  B 
and  the  objective  function  of  the  process  under  consideration  or  is  a 
consequence  of  the  general  result  that  the  optimal  policy  with  a  single 
state  variable  is  always  disjoint  1 11  • 

Bie  numerical  data  chosen  for  the  calculation  of  equations  (166),  (185)» 
and  (191)  are 


A^) 

-1 

ft* 

ho 

-  Z.Slx 

lo': 

hr-^ 

"20 

-  1.993 

xio' 

hr-1 

5,  m  10,000  Btu/lbHsol 
1^  m  20,000  Bta/lb-HBol 


n 


v^^  «  1,000  ft^/hr 

-3 
t       =  2»5  X  10  "^  hr 

I2      s   10   ft 

R      ,2—21^ 


IbHnole    R 


From  oquatlcai  (201),  the  optimal  toaperature  in  the  csm  is  calculated  as 

e^  B  930.4   R 

with  this  teniperatuare,  the  optiimua  outlet  conoaitration  ftran  the  CSTR  is 
calculated  according  to  equation  (I66)  to  be 

xj  =  0.583 

Substituting  the  above  optimal  feed  conditiai  of  the  tubular  reactor  into 
equation  (191  )•  and  solring  for  xj^^  by  trial  and  error,  one  obtains  the 
optiraal  concentration  profiles  along  the  reactor.  The  optiimun  final  laroduct 
thus  obtained  tdthout  recsytla  is 

7  »  0.916 
Ihe  optimal  taig)erature  profiles  along  the  reactor  can  be  reoorered  from 
equation  (185).  The  results  of  calculations  with  a  fractional  recycle  rate 
3  as  the  parameter  are  showi  in  Figure  14.  In  dealing  idth  equation  (191)  t 
values  of  v^  '  are  calculated  according  to  various  values  of  0  chosen. 

It  is  worth^ile  to  raantiwi  that  the  optimal  solution  obtained  can 
be  vralfied  by  the  pertuzi)ation  technique.  This  is  accomplished  by  investi- 
gatljig  the  vicinity  of  the  optimal  deoisicai  to  see  if  it  increases  the  value 
of  the  objective  functicm.  For  this  problem,  the  following  check  by 
perturbation  shown  in  Table  3  bas  been  made.  In  this  table,  the  values  of 
x^  '(L,)  in  Rows  al  and  bl  are  the  optimal  values  for  3  =  0  and  P  =  0.5 


respectively. 

It  is  also  wroth  mentiaiing  that  the  optinal  tengjerature  obtained 
Troni  equation  (201)  for  the  CSTR  is  essemtlalOy  the  saae  as  the  optinal 
tenQ)erature  at  the  inlet  of  the  tubular  reactor  without  product  recycle, 
as  given  by  equatiai  (185).  This  again  shows  the  disjoint  character  of 
optljnal  teoperature  for  the  reaction  A  —*-  B,  that  is,  the  ten^ratur* 

should  always  be  chosen  so  that  the  rate  of  reaction  will  be  as  large  as 
possible* 
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1000 


2  4  6 

Length  of  reactor^  ft 
Fig.l  4  Optimal   temperature  and   concentration  profiles   in 
a  tubular  reactor  with  product  recycle  for  the 
reaction    A^^B 

k2 


Table  3.  Perturbatiai  check  for  the  results  of  Sectioo  hd 


e^CR) 

^' 

^^^Hi^) 

aX 

930.^ 

0.583 

0.9161 

3.0 

•2 

910 

0.5815 

0.9158 

a3 

960 

o,58a;f 

0.9156 

bl 

930.4 

0.583 

0.883 

3«  0.5 

b2 

880 

0.568 

0.881 

b3 

980 

0*575 

0.882 
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NOMSHGL'.r..'RE 

A 

chad-cpj.  specie 

,(k) 

crosn  noctronal  ar'-^a  of  a  tubular  reactor  in  the  >rth  branch 

b 

order  of  inaction 

B 

cheiiical  si^ecie 

°1 

constant  in  objoctivo  functicai 

C 

chendcal  si>ecie 

^i 

concentration  of  component  i 

h 

activation  energy 

f 

number  of  final  po:^nt3 

^ 

perfomance  function 

% 

optimum  value  of  the  objective  function  for  a  N-stage  process  in 

r'ection  k-c 

g 

function  of  state  variable  associated  witli  separating  point 

% 

optirnici  value  of  the  objective  function  for  a  N-stage  process 

% 

total  amount  of  wasJiing  -rater  supplied 

h 

interval  profit  function 

H 

the  Hariiltonian  function 

k 

branch  nuriber 

h 

siJecifie  reaction  constant 

S.0 

freqttoncy  factor 

I 

independent  variable,  axial  distance  along  a  tubular  reactor 

\ 

total  length  of  the  Icth  branch  tubular  reactor 

final  value  of  independent  variable   at  the  Icth  branch 

n 

the  nth  stage,  total  nrjaber  of  branches 

H 

total  nuraber  of  stages 
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P  net  profit  function 

q  nass  flow  rate 

Q„  total  araount  of  product  H  extracted 

r  dinension  of  decision  vector  in  Section  2a  and  Section  4b, 

order  of  reaction  with  respect  to  R  in  Section  4c 

R  chenical  specie  or  z^s  constant 

R.  rate  of  production  of  ccanponeart  i 

s  dinension  of  state  vector,  order  of  reaction  vrith   resix>ct  to  R  in 

t  the  reverse  reaction  of  an  autocatalj'tic  roacticsi 

S  objective  fiuiction 

S-  objective  fiinction  of  the  ith  unit 

t  independent  variable,  tine  or  position;  holdinc  tine  in  Section 

3a  and  Section  3b 

T  final  value  of  independent  variable  t 

Tu  final  value  of  independent  variable  t  at  the  kth  branch 

u"  w"/q 

yvk;  volumetric  flow  rate  in  the  Icth  branch 

V  volune  of  a  continuous  stirred  tanlc  reactor 

^  irasliinc  crater  flow  raoo  into  the  xi^'a   stage 

X  state  vector 

Xj,  no?-ar  conc-^ntntion  of  ccxiponent  i  in  '^oct-'.on  3a  Jind  Hoction  3b 

V**  composition  of  errtmct  at  oi.itlot  of  the  nth  extractor 

y  optinan  value  of  a  state  variable  at  the  end  of  a  tubular  reactor 

2  adjoint  vector 


Greek  Letters 

tt    initial  value  of  a  state  vector  in  Section  ^;  proportionality 

constant  of  equilibrium  relation  in  Section  'tc;  fracticmal  amount 
of  feed  stream  in  Section  4d 

«f    Initial  condition  of  state  variable  x 

i.  i 

0  fractional  amount  of  recycle  stream  in  Section  4d 

P-  initial  condition  of  x^  in  the  intermediate  feed 

Y  volumetric  flow  rate  fractic«  of  the  feed  stream 

Y  volumetric  flow  rate  fraction  of  the  intenaediate  feed  stream 


I 


dimensionless  holding  time  in  the  ith  reactor 


vk)  (k) 

Tr  closer  approxlraaticn  to  the  correct  values  of  x£  ' 

9  decision  vector 

f  the  optimum  value  of  6 

5^^  better  approximate  rale  of  i^H\)   In  the  Runge-Kutta  method 

p^  density  of  a  mixture 
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The  main  purpose  of  this  vork  is  to  s^ow  how  the  maxixnm  principle  and 
djmaraic  prc^rannning  can  bo  applied  to  topolqjically  complex,  chaoical  rate 
processes,  i«urticularly  to  flow  type  oh«aioal  reactor  systems.  The  study 
i^ich  treats  four  escanipiLes  in  detail  consists  of  t^  parts* 

The  first  part  is  an  investigation  of  optimal  tec^perature  profiles  in  a 
series  of  plug  flow  reactors  with  either  fresh  feed  fed  to  or  intermediate 
product  withdrawal  from  an  internal  point  by  using  the  laaxiHua  principle. 
For  both  exajsples*  it  is  assuned  that  a  first-order  consecutive  reaction 
A — »B — ^C  is  taking  place  in  the  reactors.  The  split  type  two-point- 
boundary  value  problems^  tWich  are  encountei^,  are  solved  numerically  by 
using  an  iterative  Runge-Kutta  inethod  and  backward  int ^ration.  In  the  case 
of  intenaediate  product  withdrawal,  it  is  found  that  the  optimal  temperature 
profile  in  the  second  reactor  is  almost  flat,  i.e.  isothermal  operatiao  in 
the  second  reactor  will  meet  the  optimum  ccsidition.  For  the  case  of  inte)> 
nediate  fresh  feed,  the  effects  of  fractional  amounts  of  the  intermediate 
fjpesh  feed  on  the  optimal  temperature  profile  are  also  studied. 

The  seccaid  part  deals  with  the  combined  use  of  the  raaxlmam  principle 
and  dynamic  programming.  General  working  equaticas  and  conixxtatianal  schemes 
for  a  sequ«itial  union  of  these  two  raethods  are  developed.  It  is  found  that 
the  coBfcined  use  of  these  two  modem  optiniization  techniques  offers  many 
advantages  over  either  caie  of  these  two  methods,  particularly  to  the  geomet- 
rically cooplaoc  heterc^eneous  processes.  To  illustrate  the  use  of  these 
equations,  two  examples  are  woi^ced  out  in  detail* 

m  the  first  example  a  tubular  flow  reactor  followed  by  a  multistage 
cross-current  extraction  process  is  considered.  An  autocataljrtic  waction 
is  considered  in  the  reactor.  The  working  equations  are  derived  for  an 
objective  fuiwtion  i^ich  maximiaes  total  net  profit.  This  profit  is  maximized 


1^  the  optiraal  choice  of  the  terapearat'.ipe  profile  in  the  tubular  reactor  and 
hy  the  optirsal  allocation  of  washing  water  in  the  cross-current  orbraction 
towers* 

The  second  exaraple  consists  of  a  continuous  flow  stirred  tank  reactor 
followed  by  a  tubular  flow  reactor  with  product  recycle  to  an  internal  point* 
A  first  order  reversible  reaction  is  assuaed  to  take  place  in  the  reactors. 
Optimal  teapejratures  in  the  Csra  and  along  the  tubular  reactor  are  determined 
in  order  to  raaxLiniae  the  concentration  of  the  desired  final  product,  15ie 
effects  of  recycle  to  an  intenml  point  are  also  investigated. 


